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ABSTRACT 



One of the responsibilities of the personnel planner of the Bureau 
of Naval Personnel is the projection of current personnel inventories 
into the future and the determination of the personnel replacements 
that are necessary in order for these inventories to meet stated re- 
quirements. The personnel planner presently has a method of 
solution to this problem no more sophisticated than hand computations, 
and consequently it is a long and arduous task. 

The process by which a personnel inventory changes with the 
passage of time is suitable for programming on a digital computer. 

The complexities that this problem generates because it considers 
the behavior of people only serves to amplify the utility of a digital 
computer to perform this function in preference to hand computations. 

This thesis is the development of the program PIPE, written in 
FORTRAN 60 for the 1604 digital computer, which will project current 
enlisted personnel inventories of certain subsystems of the Navy for 
as much as five years into the future, and determine the least 
expensive personnel replacements necessary to enable these pro- 
jected inventories to meet stated requirements. The program reduces 
the solution time of this particular problem of the personnel planner 
to minutes. 

PIPE is not designed to be a mechanical process to be used 
blindly. Rather, this program can be used effectively to improve 
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the art of personnel inventory projection when the personnel planni 
uses his experience to properly interpret the results. 

An analysis of PIPE with respect to its current and future 
potential to the personnel planner led to certain conclusions and 
recommendations that are enumerated in this paper. 



PREFACE 



The Naval personnel planner is faced today with an increasingly 
complex problem: how to provide enough adequately trained enlisted 
men, on time, to meet the requirements of rapidly expanding new 
weapons systems. 

The authors became interested in this problem during the 
course of temporary additional assignments in Washington, D. C. , 
during the summer of 1963. These assignments were made as part 
of the Operations Analysis Curriculum at the U. S. Naval Post- 
graduate School. Lieutenant Commanders Meeks and Metcalf spent 
six weeks at the Bureau of Naval Personnel, and during this period 
became well acquainted with the Bureau's current methods used in 
meeting this problem, as well as with some new concepts that the 
Bureau was considering for use in this problem area. Lieutenants 
Greenhalgh and Harper were working at the Special Projects Office, 
and there became impressed with the concern shown for the dif- 
ficulties involved in planning for sufficient enlisted personnel to 
meet the requirements of the Polaris program. 

The four officer -students thought it to be beneficial to pool 
their experience and knowledge of the problem in a combined effort 
which would attempt to derive a simple yet useful program for the 
solution of this personnel problem. Upon closer examination it was 
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seen that the problem was large and complex enough to indeed merit 
the combined efforts of the four officers, and that such effort, if 
carried to a conclusion comensurate with their proposals, would 
satisfy the requirements for a thesis in the Master's Degree Program 
of the Operations Analysis Curriculum of the Naval Postgraduate 
School, 

The result of this effort is Personnel Inventory Projection, 
Enlisted, (PIPE). PIPE consists of a linear programming model 
along with a digital computer program that will project personnel 
inventories and determine the least expensive personnel replace- 
ments that are necessary to enable these projected inventories to 
meet stated requirements. 

The development of the model and the computer program has 
been subdivided into sections so as to allow readers with different 
interests to read only the pertinent sections of the thesis. 

A detailed discussion of the problem of projecting personnel 
inventories and the reasons why the authors were motivated to 
undertake this problem can be found in Chapter I, 

The derivation of the mathematical equations which serve both 
to project a current inventory and to provide the necessary infor- 
mation for using a linear programming technique for computing 
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least expensive personnel replacements, can be found in Chapters 
II and III, and Appendix B. 

The assumptions considered necessary in the derivation of 
the mathematical equations can be found in Appendix A. 

The procedures utilized in writing a digital computer program 
for the mathematical equations and the linear programming technique 
can be found in Chapter IV. 

Instructions for the use of the computer program with a de- 
tailed description of the necessary inputs and the resultant outputs 
can be found in Appendix D. A program listing can be found in 
Appendix E. 

The program was run with a representative current inventory 
with various values of the input parameters, and a sensitivity analysis 
conducted upon the results, A summary of this analysis can be found 
in Appendix C. These results were obtained using the U. S. Naval 
Postgraduate School's CDC 1604 computer. 

Conclusions and recommendations can be found in Chapter I, 
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CHAPTER I 



NON-TECHNICAL DISCUSSION 
A. Introduction 

I. Definition of the Problem 

As the modern Navy becomes more technically oriented, 
the personnel planner is faced with an increasingly complex problem: 
how to provide adequate trained enlisted personnel, on time, to 
meet the requirements of rapidly expanding new weapons systems. 
While similar problems have always confronted the planner in 
varying degrees, today's management decisions concerning per- 
sonnel are more complex than ever before. 

The Navy has a serious shortage of qualified manpower 
for its technical programs. In the past, with the criticality of talent 
not so acute, the planner could afford some margin for error in his 
planning. But today he must optimize the utilization of the scarce 
resources under his control, or soon find a reduced effectiveness 
of these critical programs due to a lack of qualified personnel. 

The technical programs require long lead times for 
training. This further complicates training planning because input 
quotas must be determined further in advance, school sttrition rates 
must be more carefully analyzed, and more accurate accounting of 
personnel in training must be attained. 
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The high turnover of personnel in critical skills requires 
frequent and increasingly costly replacement. A large part of each 
man's first enlistment in a critical rating is spent in school or some 
other form of training, and reenlistment rates for these personnel 
are discouragingly low. To achieve maximum return on this high 
training investment, there must be a most careful analysis to 
determine the optimum employment of each man in his particular 
skill. 

Any technical personnel training plan should take into 
account the future needs of the whole Navy, since the current 
technical programs may not always have the high priority that they 
now enjoy. Thus, any personnel planning must remain flexible and 
carefully consider the future needs of the Navy as they become evi- 
dent. 

For flexibility and reliability in planning, a fast and 
efficient method of changing the countless variables involved in 
personnel management is required. The personnel planner, using 
hand computations alone, is able to consider only one or two alter- 
natives; he then implements the selected plan, and relies on later 
revisions to make the program suit the Navy's needs. The result 
is a loss of time, talent and efficiency. 

The choice of any training plan should be based not only 
on its output of trained personnel, but also on its side effects on 
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such factors as morale, reenlistment rate, cost and utilization of 



facilities. Only in the optimization of all these factors can the 
planner employ his personnel resources with maximum efficiency. 

The general approach taken in this paper is to construct 
a mathematical model of a "critical" personnel subsystem, and then 
to program it on a digital computer. The objective is to provide the 
personnel planner with a "tool" with which he is able to predict 
future inventories from a given starting inventory, under known re- 
quirements and certain acturial factors. In the process of obtaining 
the desired projected inventories of the initial group, optimal 
training inputs are determined. 

2. Background 

The Operations Analysis Curriculum of the U. S. Naval 
Postgraduate School requires that each student write a thesis during 
the second year of study. Two such theses written in 1963 were 
instrumental in stimulating the authors' interest in the broad area 
of personnel problems. They were Operations Research and 
Management of the Navy's Personnel System by LCDR John W. 
Walden, USN, j^lj ,*and An Analysis of the Cost and Requirements 
of the Fleet Ballistic Missile Submarine Personnel Subsystem by 
LT Edwin M. Baldwin, USN, [z] . 

^Numbers in brackets refer to references in the Bibliography, 

Page 89- 
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LCDR Walden's thesis essentially was a broad look at 

how operations research techniques might be applied to the Navy's 

complex personnel system, and cited various operations research 

type studies concerned directly with military personnel management. 

One of LCDR Walden's objectives was to stimulate further research 

in the area of applying OR techniques to personnel problems. Among 

his recommendations, he stated: 

a necessary first step has to do with studying the 
(personnel) system in its entirety. Some model of 
our personnel system must be developed and some 
meaningful measures of effectiveness derived. 

The authors were further motivated to investigate 
the Navy's personnel problems by their summer field trips to 
Washington, D. C. , during the period 1 June to 12 July 1963, as 
part of the Operations Analysis Curriculum. LCDR's Meeks and 
Metcalf were assigned to the Bureau of Naval Personnel, while 
LT's Greenhalgh and Harper were attached to the Special Projects 
Office. In the course of the summer assignments, the officers 
assigned to BuPers became aware of the compelling need for the 
Navy's personnel planners to be able to project personnel inven- 
tories under the constraints of varying policy changes and special 
programs. The officers assigned to the Special Projects Office 
gained an appreciation for the increasingly acute personnel problems 
associated with the rapidly expanding demand for the critically short. 
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highly qualified, highly trained personnel of our FBM program. 

While in Washington, the authors attended a briefing 
at the Special Projects Office where a paper entitled Nuclear Officer 
Training Requirements Program by CDR Charles E. Woods, USN, ^3 
was presented. Briefly, this paper was concerned with the opti- 
mization of the numbers of nuclear -trained officers in various 
seniority groups to meet the demands of the expanding SSBN and SSN 
building programs. 

This study was of paramount importance in the recent 
reevaluation of the training plans for officers entering the SSBN 
and SSN programs. CDR Woods, a graduate of the Operations 
Analysis Curriculum of the U. S. Naval Postgraduate School, 
utilized many operations analysis techniques in his study. His 
paper was therefore of interest because it represented an example 
of how these techniques could be successfully applied to the officer 
personnel system. It was also of interest for the possible appli- 
cation of a similar model to the enlisted personnel system. 

Out of a recognition of the vital importance of the 
personnel subsystem to any weapons system, a desire to apply 
operations research techniques to the Navy's personnel system, 
and a mutual interest in "people", the four authors were motivated 
to undertake this paper. 
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Concept of PIPE (Personnel Inventory Projection, Enlisted) 



1. Objective of the Model 

The Navy personnel planner needs a method which will 
quickly, accurately and with broad flexibility: 

a. Project inventories of personnel for a desired time 



period; 

b. Compute optimal training inputs under known re- 
quirements and current acturial factors; and, 

c. Predict excesses and shortages of the projected 
inventory by pay grade. 

The objectives of PIPE are to accomplish the above with 
a digital computer program which will require a minimum of running 
time, be simple to operate and modify according to the user's re- 
quirements, and employ imputs in a form readily available to the 
user. In other words, it is hoped that, by using this model, the 
personnel planner can get a rapid and reasonable estimate of the 
personnel situation of the future, and experiment with ways of 
improving this situation as he sees fit. 

2. Classification of Naval Enlisted Personnel 

The structure of the Naval enlisted personnel system is 
complex and varied, with great differences between men as to 
experience level, skills, and criticality. In order to predict the 
number of men in this structure at some future date, it is first 
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necessary to examine the personnel system for some kind of con- 
sistency of qualities so that the men can be grouped and examined 
with respect to recurring phenomena. A breakdown by ratings alone 
may not be sufficient for this purpose, because some ratings may 
include a broad spectrum of specific skills, some of which may be 
in critical supply, and some of which may not. 

Within the critical ratings, however, the Naval Enlisted 
Classification Code (NEC) does provide such a grouping of personnel. 
Men with a common NEC share comparable skills and criticality; 
they also share approximately the same attrition, reenlistment, 
retirement, and promotion factors. Therefore, the enlisted men 
considered in this model may be delineated by NEC in the critical 
ratings, or by any other criteria that provides similar character- 
istics. 

3. The Inventory Matrix 

The men in a specific group (i. e. , with common acturial 
factors) are next further subdivided by both pay grade and years of 
obligated active service. This delineation lends itself to a matrix 
formulation of the initial inventory. The matrix consists of cells 
representing all personnel in the group, with the numbers in 
appropriate cells representing personnel according to their present 
obligated service and pay grade. Within each pay grade, the matrix 
is further categorized by specified enlistments. This allows 
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personnel in the inventory to be delineated by total years of active 
service as well as by pay grade and obligated service. In the model, 
total years of service are analogous to current enlistment, where it 
is assumed that all enlistments are of six years duration. 

Hereafter, the word "inventory" will always denote this 
type of matrix description of the particular group of personnel under 
consideration, with cells of the matrix representing numbers of 
personnel according to their obligated service, pay grade, and 
current enlistment. This type of formulation is considered necessary 
because this model deals with people , and as such, a means must be 
provided for accounting for the differences between men within the 
group. It is felt that such a matrix expresses the differences in 
professional capabilities of the personnel in the group. 

4. Change of Inventory With Time 

As indicated above, a given cell of the inventory matrix 
represents the number of personnel categorized by: 

a. Obligated service 

b. Pay grade 

c. Current enlistment 

With the passage of an increment of time there will be 
a change in the inventory matrix because at least one of the categories 
above will change for all personnel. Thus, a man who is represented 
in the inventory of a given cell at a certain time will necessarily be 



8 




i 



represented in a different cell after an increment of time. This 
change of the inventory with time, then, can be visualized as a 
"movement" of personnel from one cell to another within the in- 
ventory matrix. Throughout the remainder of this paper, the 
"movement" of personnel within the inventory matrix will refer to 
this concept of the change of the inventory with time. 

Because personnel in the critical ratings require 
approximately one year of formal schooling or other training before 
they are ready to join the operating forces, it was decided to look at 
time in discrete increments of one year. 

There are four basic acturial factors or rates, which 
are a characteristic of any particular inventory in the personnel 
system. They are defined as follows; 

a. Attrition rate: a percentage of men who leave the 

inventory due’ to death, sickness, disciplinary action, lack of 
aptitude, efc. 

b. Retirement rate: the percentage of men who leave 
the service by retiring after 20 years or more of active service. 

c. Promotion rate: the percentage of men who are 
promoted from any one pay grade, to the next higher pay grade. 

d. Reenlistment rate: the percentage of men who 
reenlist in the Navy for a period of six years and thus remain in the 
inventory. 
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It is the application of these actuarial factors in the approp- 
riate manner which causes the change in the basic inventory per unit 
time. Thus, it can be said that these four rates cause personnel to 
"move" through the inventory matrix. In the sense that these factors 
cause a change in the inventory matrix with time, they can be 
thought of as "operating" on the basic inventory matrix to "transform" 
it into a new, or changed inventory. Accordingly, these basic actu- 
arial rates of attrition, reenlistment, retirement, and promotion 
will be referred to as "operators". 

Upon examining these four operators as they affect the 
inventory matrix, the follpwing conclusions result: 

a. The attrition operator will affect every cell in the 
inventory matrix, serving to reduce the number of people in each 
cell as time passes. 

b. The retirement operator will affect only those cells 
which represent personnel with 20 years of active service or more. 

c. The promotion operator will apply to all cells 
which are determined to contain men eligible for promotion. 

d. The reenlistment operator will apply only to those 
cells which represent men in their last year of obligated service. 

It is readily seen that the "operators" must be applied 
to different cells in different ways in order to cause movement 
through the inventory matrix in a realistic manner. Thus each 
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operator takes on the form of the inventory matrix, whereby each 
cell of the inventory matrix has a corresponding attrition, retire- 
ment, promotion, and reenlistment rate associated with the respective 
operator matrices. Through a proper means of indexing, the 
appropriate operators are thereby applied to their corresponding 
inventory cells, and the personnel represented in the inventory matrix 
will be "moved" through the matrix with the passage of time. 

Due to the action of the operators, the numbers in all 
cells of the inventory matrix will be reduced to reflect attrition, 
and the numbers in certain cells will be reduced to reflect retire- 
ment and non-reenlistment. To replace these men there must be 
a means for inserting new men into the program in any given year. 

This is done with training inputs, which again take the form of a 
matrix corresponding to the inventory matrix. In this way numbers 
in selected cells of the inventory matrix can be increased and these 
cells can be varied to reflect the user's choice of training inputs. 

5. The Projection Equations 

The passage of time: 

a. Generates movement of personnel through the cells 
of the inventory matrix under the influence of the operator matrices; 

b. Provides for the addition of personnel as training 

inputs . 
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With this knowledge, projection equations can be derived 



that express the projected inventory in terms of the original inventory, 
operators, and training inputs. The only unknowns of these equations 
are the training inputs. 

In order to utilize these equations, personnel require- 
ments for the future must be made available. Basically, these 
requirements emanate from the Department of Defense in The Five - 
Year Force Structure and Financial Program . This model uses 
these requirements grouped by pay grades for each year of projection. 

If the requirements are set as linear bounds on the pro- 
jection equations, then by a mathematical technique known as linear 
programming, the unknown training inputs can be determined in 
such a way as to satisfy the requirement in a manner which is 
optimal with respect to cost and/or numbers of personnel. 

With the solutions to the optimal training inputs avail- 
able, the projection equations then can be used to express explicitly 
the projected inventory. 
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C . Other Personnel Computer Models 

Three separate computer models generally concerned with the 
projection of personnel were studied prior to the conception of the 
models developed in this thesis. A brief description of these systems 
is considered appropriate in order to place the subject model in 
per spective. 

1. Project MOON 

Project MOON (Meeting of Operational Needs) is an 
enlisted personnel simulation system that was developed by the 
American Institute of Research under contract from the Personnel 
Research Division (Pers 15) of the Bureau of Naval Personnel. 

MOON provides for the projection of a selected rating for any period 
of from one to five years into the future. These projections are 
made either (1) on the basis of the continuance of current policies, 
planning methods, and actuarial factors, or (2) with the introduction 
of selected changes in any of these policies, planning methods, or 
actuarial factors. Thus, MOON can be used to estimate the future 
state of a selected rating under current conditions, or test the 
effects of the changes in these conditions over the desired period of 
projection. 

One of the essential features of the MOON program is 
that it operates on "real" people; the basic starting sample popu- 
lation consists of extracted service records taken from the BuPers 
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magnetic tape master file. Each "person" in the sample population 
is denoted by a sequence of alpha-numeric characters that denote 
nine basic simulated personnel variables which are: 

a. Pay grade 

b. Primary Naval Enlisted Classification Code (NEC) 

c. Secondary NEC 

d. ETST Score (or MECH) 

e. USN/USNR 

f. Time in grade 

g. Obligated active service 

h. Length of service 

i. Sum of the basic test battery scores, (GCT, ARI, 
and ETST or MECH) 

Other inputs to the program consist of fourteen directly 
controllable variables which represent the various policy and plan- 
ning methods previously mentioned, and seven indirectly controlled 
variables which represent the actuarial factors involved. 

The main program of Project MOON simulates the 
passage of time by performing various peripheral routines on the 
appropriate sample population. The actual people, (i. e. , extracted 
records), considered by the main simulation program are selected 
randomly by means of Monte Carlo techniques from the population 
sample. 
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The basic input of each projection run is a set of tabu- 
lations representing the state of selected ratings at the start and at 
the end of each fiscal year of the projection period. These reports 
are cross -tabulated in various ways, producing an awesome array 
of data, including various histograms. 

In summary, Project MOON is a highly sophisticated, 
complex program which is considered to have two particular limi- 
tations: 

a. The running rime (computer time) per year of 
projection for one selected rating is about one hour; hence, a 
five-year projection requires about five hours of computer time. 

This is a function of the complexity of the program, particularly 
with regard to the large number of random selections made in the 
program. 

b. Both the annual updating procedure and the annual 
validation procedure for MOON require massive manual data- 
handling efforts. 

2. The Nuclear Officer Training Requirements Program 

The Nuclear Officer Training Requirements Program ^3 
is the computer model which was developed from the mathematical 
model conceived by CDR Charles E. Woods, USN, who was then 
attached to the Office of the Chief of Naval Operations (Op-31). 
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Commander Wood's model was programmed by Booz -Allen Applied 
Research, Inc. , Washington, D. C. , under contract from the Special 
Projects Office. 

This model uses elementary matrix methods to predict, 
year by year, the number of nuclear trained officers in the various 
seniority groups, based on assumed inputs into the training program 
and estimates of attrition for the different seniority groups. By 
varying the input parameters, the effects on the distribution of 
officers in various seniority groups can be examined. The predicted 
numbers by seniority groups are then compared with known require- 
ments. Estimated requirements for trained officers are based on 
known shipbuilding schedules and established or postulated seniority 
criteria. The inputs of officers entering training for successive 
years are treated as unknowns. 

This model uses linear programming techniques to 
minimize the total number of officers requiring nuclear training 
based on the requirements, seniority criteria and attrition factors. 
The outputs of the computer model is training requirements by 
seniority groups per fiscal year. 

3. SMS Personnel Projection Model 

The Surface Missile Systems Personnel Projection 
Model is an elementary, deterministic computer model which was 
developed for the Special Surface Missile Systems Task Forces 
(Pers 17) by the computer facility at the Naval Weapons Laboratory, 
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Dahlgren, Virginia. Although the authors were able to obtain a 
program listing for this model, the model is undocumented as far 
as basic explanation and concept are concerned. 

Prior to the development of the computer model, Pers 17 
undertook the major task of purifying its basic inventory of trained 
personnel by attempting to identify adequately all personnel who had 
ever been trained for one of the critical NEC's involved in the sur- 
face missile systems. With their inventory purification complete, 
they next developed a simplified model which would project their 
inventory into the future under certain known acturial rates such 
as attrition, reenlistment and retirement. The program's output 
consists of predicted totals of personnel in various NEC's for each 
six -month increment into the future. 

D. Comparison of Models 

PIPE is considerably less sophisticated than Project 
MOON; it does not consider as many variables, and it does not 
attempt to project the inventory of an entire rating. It does, 
however, consider the important factors in the projection of per- 
sonnel inventories, and because of its simplicity it requires only 
a few minutes to run on the computer. Also, it does not require 
the massive accumulation of data that MOON requires. 

The SMS projection model of Pers 17 treats an NEC as 
a single entity, but does not consider requirements or training inputs. 
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PIPE is more sophisticated in that the NEC's are delineated by pay 
grade and are compared to requirements to obtain optimal training 
inputs for every year of projection. 

The Nuclear Officer Training Program model was de- 
signed to project inventories of officers; it cannot be directly trans- 
ferred for use as a projection of enlisted men due to the added com- 
plexity of the enlisted personnel system. But just as CDR Wood's 
model projected inventories of officers by submarine billets over 
a period of time, so PIPE will project inventories of enlisted 
personnel by pay grades in NEC's or similar groups. Thus, there 
is a similarity between the two models in concept and approach, 
but PIPE is considered more complex due to the greater complexity 
of the enlisted personnel system. 
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E. 



Conclusions 



1. Linear programming is a technique traditionally associated 
with the optimal allocation of resources with respect to minimum 
cost and/or other requirements. PIPE demonstrates that linear 
programming can provide valuable assistance to the personnel 
planner concerned with the projection of current inventories into 
the future, the necessary personnel inputs over the years of pro- 
jection, and the problem of ensuring that the projected inventories 
meet the requirements of the future. PIPE specifically provides 
estimates of the annual training inputs necessary to meet the future 
requirements of critical rating/NEC groups at minimum cost, where 
costs are measured either in number of personnel or monetary units. 

2. PIPE is a unique program for the projection of personnel 
inventories because it considers requirements five years in advance 
of the present time in deriving its minimum cost solution. There- 
fore, PIPE is a long-range planning tool in that it "looks ahead" to 
the future and computes the best solution to the problem considering 
the requirements for personnel of each year for five years. 

3. The most difficult part of any personnel inventory projection 
of an enlisted personnel subsystem is the initial breakdown and 
categorizing of the personnel in the subsystem. While generali- 
zation leads to inaccuracies, it is also true that great detail brings 
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about intricately complicated and time-consuming computer programs. 
PIPE is sufficiently detailed to be useful to the personnel planner, 
yet basic enough to require simple inputs and a minimum of com- 
puter time. 

4. It is recognized that many complex prediction equations have 
been developed in this thesis, and that no matter how logically they 
have been derived, they can not be considered dependable until their 
solutions are exposed to the realities of empirical data, PIPE was 
operated with actual input data from BuPers for Fiscal Year 1963; 
it cannot be validated until historical data becomes available. 

5:. The use of obligated service as a criterion for the accounting 
of personnel in the inventory allows an accurate and efficient use of 
reenlistment rates as input parameters in the projection of personnel 
inventories. No matter what the length of a man's time in service, 
or his pay grade, or his enlistment, he must either reenlist or 
leave the service during the last year of his enlistment. The obligated 
service criterion pinpoints him in the inventory in such a manner that 
the reenlistment rate can be directly applied to him when he is in 
his last year of obligated service. PIPE reflects the use of this 
criterion. 
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6. The obligated service criterion allows training input to be 
delineated further than has heretofore been the case. Training 
inputs to the subsystem can therefore be specified by their length 

of service, pay grade, and obligated service, thereby allowing more 
specific control over these inputs to a personnel projection. 

7. Experience has shown the personnel planner that controlled 
variations of promotion, retirement, reenlistment, and attrition 
rates cause predictable trends in the training inputs and associated 
costs necessary to meet future personnel requirements. Using the 
actual current inventory of a critical rating, with various values 

of the input parameters, (promotion, retirement, reenlistment and 
attrition rates), the resultant outputs of the model were analyzed 
with respect to the sensitivity of the solutions to the actual values 
of the input parameters. The details of this sensitivity analysis 
are described in Appendix C. This analysis supports the afore- 
mentioned trends, and hence the model is considered to use these 
rates effectively in the projection of an inventory; it shows that the 
model reflects for all cases considered the proper direction of the 
trends anticipated for the various input parameters. 
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F. 



Recommendations 



1. Test the validity of the model using historical data. When 

the information becomes available for personnel inventories separated 
by at least one year, then the model should be tested for its approxi- 
mation of the real life situation. 

2. Investigate the possibility of computing reenlistment, retire- 
ment, and attrition rates from information on the BuPers master 
magnetic tape files. 

3. Investigate the actual distribution of training input by pay grade, 
obligated service, and length of service, in order to evaluate the 
validity of the uniform -distribution assumption used in the model. 

4. Operate the model with various training-input distributions, 
and analyze their effects on the solution with respect to the total 
training input, and excesses and shortages in the different pay grades. 

5. Investigate the various methods by which the relative costs 

of the training inputs from the different pay grades can be computed. 
Since these costs serve to control the linear programming solution, 
the accuracy with which these costs are determined will directly 
affect the utility of the model. 
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6. The linear program in this model uses restraint equations 
with many lower limits, while the only upper limit represents the 
permissible excess of inventory over requirements. Investigate 
the feasibility of placing additional upper limits on the restraint 
equations. These upper limits could correspond to the influences 
of such factors as the capacity of training facilities, and the availa- 
bility of personnel in certain pay grades. 

7. With a very minor program change, which would enable per- 
fectly arbitrary training plans to be an input, project current inven- 
tories with a proposed training plan and compare the projected 
inventories with requirements. This procedure enables the person- 
nel planner to compare numerous training plan proposals. 
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CHAPTER II 



BASIC MODEL 



A. Concept 

The basic problem being undertaken is one of complex inventory control, 
characterized by the following facts: 

1. Items in the inventory are people. 

2. People with different professional capabilities must be subdivided to 
provide distinction. 

3. The number of people to be foimd in the inventory is affected by numer- 
ous factors, the values of which depend on the particular people imder 
consideration. 

4. During the passage of time, the people assume new potentials, and this 
change can occur in various ways. 

The initial problem is that of creating a method of tabulating the inventory in 
a manner which is compatible with the characteristics enumerated above. The 
method utilized by the Officer Model discussed in Chapter I uses "years of com- 
pleted service" to subdivide the inventory as to professional capabilities. This 
characterization is inadequate for enlisted men. The criterion for selection of 
subdivision rules must be such as to provide the user of the model with the maxi- 
mum usable information concerning the professional potential of the inventory. 

The term professional potential as it applies to the enlisted personnel inven- 
tory requires definition. The various billets required to be filled by a given inven- 
tory vary in complexity and demand different levels of education and experience. 
Additionally it is important that an inventory specify the expected tenure of service 
of each man. Professional potential is therefore defined as a function of educa- 
tional level, experience level, and expected tenure of service. It is obvious that the 
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single parameter of "years of completed service" as used in the Officer Model is 
inadequate to define sufficiently the professional potential of an enlisted man. 

Cdr. E. Napier, of the Navy's Special Projects Office, suggested subdividing 
the inventory by enlisted pay grade and years of obligated service. Such a repre- 
sentation is shown in Figure II-l. This subdivision, although satisfactory at first 
glance, does not provide for sufficient breakdown by experience level. This can 
best be illustrated with an example. 
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Figure II-l 

Subdivision of Inventory 



Consider 100 men of pay grade E-6 with three years obligated service. It 
would be far more meaningful to state that of these 100 men there are 15 men with 
about two to four years service, 60 men with about ei^t to ten years service, and 
25 men with about 14 to 16 years service. This information provides a much more 
informative description of the experience level of the inventory and allows the 
personnel administrator to appreciate more fully the professional potential of the 
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total number of men in any one pay grade. 

This then leads to the subdivision of the inventory matrix as indicated in 
Figure II-2, wherein the columns within any one pay grade are designated as 
successive enlistments. It will be noted that the number of columns, and hence 
enlistments, which are associated with each pay grade varies. By analyzing repre- 
sentative current inventories of several critical ratings through information pro- 
vided by the Bureau of Naval Personnel, it was determined that personnel on 
certain enlistments are not expected to be found in certain pay grades. This infor- 
mation dictated the specific matrix as shown in Figure II-2. 

In order to use this inventory array it is necessary to create some notational 
scheme to represent each cell in the array. If this problem is approached by 
treating the array as a conventional matrix and using double subscripts to identify 
each cell within the array, the result proves to be overwhelming in notational 
complexity. This leads then to the scheme of assigning a single subscript to each 
cell in the order indicated in Figure II-3, which simplifies the notational difficulties 
immeasurably, and more importantly, provides greater flexibility in the ultimate 
writing of a digital computer program for the mathematical model. This array is 
described as a 6x12 matrix with 72 cells. 

With the passage of time, the number of people to be found in any cell will 
change providing there is no additional input. The reasons for this change can be 
attributed to four factors; 

1. attrition rate 

2. retirement rate 

3. promotion rate 

4. reenlistment rate 

In the model, these factors will be represented by so-called matrix operators. 
The subdivision of the inventory array by cells demands that a separate set of 
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Figure II-2 
Inventory Matrix 
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Figure II-3 

Subscripted Inventory Matrix 
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values for these operators be provided for each cell. The following definitions are 
therefore necessary: 

n = the number of the cell in the matrix under consideration 
(n E 1, 2, 72) 

a = attrition rate for cell number n 
n 

OL - [a^] = complete 6x12 array of attrition rates 

8 = retirement rate for cell number n 

n 

jS = “ complete 6x12 array of retirement rates 

y ^ = promotion rate for cell number n 

y = ^ complete 6x12 array of promotion rates 

6 = reenlistment rate for cell number n 

n 

6 = [6^] = complete 6x12 array of reenlistment rates 

By these definitions the necessary specification of any operator demands an array 
similar to the inventory matrix. For example 



a = * 



With these definitions established it is possible to consider the ensuing move- 
ment of people from one cell to another during the passage of time. In order to 
determine how to represent the change in inventory during a passage of time it is 
helpful to specify the conditions imder which a change in inventory can occur. These 
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conditions are threefold and are specifically: 

1. The passage of one year permits a man to move to the next cell in the 
array, subject to the restrictions of attrition, retirement, and reenlist- 
ment; 

2. The movement described in (1) plus an additional move into the next 
higher pay grade by promotion; 

3. The addition to a cell from outside the array as a result of a training 
input. 

The nature of the inventory array dictates clearly which of conditions (1) or (2) 
can prevail for each and every cell. These are illustrated in Figure n-4. There 
exist two basically different flows, the existence of which can be appreciated from 
Figure n-4 as one notes that the line "due to promotion" is discontinuous; the first 
part of which is associated with cells n*8 to 13, the second part of which is associ- 
ated with cells n*20 to 61. Specifically this discontinuity is due to the varying 
number of enlistments to be found in any one pay grade, and hence the varying 
number of cells to be included in a jump from any one pay grade to the next. The 
cells for which there is no corresponding point on the line "due to promotion" are 
therefore cells to which one cannot be promoted, an example of which is n*16. The 
existence of these cells can be appreciated by observing exactly what enlistments the 
cells are part of and then the non existence of the enlistments from which their inputs 
would have had to come. There also exist in Figure II-4, cells for which there can 
be no input except a training input, an example being n=7. Every cell can receive a 
training input, condition (3), and is therefore not included in Figure II -4. 

These concepts necessitate the following definitions: 

X = original inventory matrix (6x12) 

= original inventory in cell number n (n* 1, 2, . . . , 72) 
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CELL NUMBER AFTER ONE YEAR 
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P(m) 




projected inventory matrix (6x12) at the end of the m^^ year 
of projection 


Pn(m) 


= 


projected inventory in cell number n at the end of the m^ 
year of projection (n» 1, 2 , . . . , 72) 


P(m) 


= 


[ Pn(m)] 


T(m) 




training input matrix (6x12) during the m^ year of projection 


T^(m) 


s 


training input into cell number n during the m^^ year of 
projection (n=*l, 2, ..., 72) 


T(m) 




[ T^(m)] 



The equations which will be written in terms of these quantities will ultimately 
be programmed for a digital computer. Since these programs must be written in 
capital letters, the symbols used above to define the cells of the matrices are capital 
letters vice the more common usage of lower cased letters. 

With these definitions established the next step is to construct the mathematical 
equations that will describe the projected inventory matrix, cell by cell (Pj^(m)), as 
a function of the cells of the original inventory matrix (X^), the cells of the operators 
(q^, 6jj), and the training input for each of the m years of projection (T^(m)). 

The initial time period to be considered is one year (m* 1). 

When these equations have been derived, a comparison will then be possible 
between the projected inventory and the stated requirements. If the training inputs 
for each cell (T^(l)) are unknowns, then the projection equations, when compared 
with the requirements, yield the necessary inputs to a linear program, a mathe- 
matical technique that will minimize the training inputs that are necessary such that 
the projected inventory (P(l)) meets the requirements. Having found the solution for 
the only unknowns in the projection equations, (training inputs), the projected inven- 
tory can then be e^ressed explicitly. 

B. Formulation 

The equation representing the inventory of each cell in the projected inventory 
matrix depends on which cell is being considered. Therefore, it is helpful to repre- 
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sent the components of the equation for each cell with a flow chart. 
Figure II-5 shows the four possible inputs to any cell, Pn(l)» of 
the projected inventory after one year of projection. For a given 
cell, Pn(l)» only a subset of these inputs are applicable. For 
this reason each subsection of Figure II-5 is explained in detail: 




Projection Flow Chart 
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Figure II-5A 

If promotion is possible from cell Xjj_i 3 of the original 
inventory, then the inventory in cell X^-lS must be subjected 
to attrition and retirement rates. That part of Xjj _]^3 that survives 
these rates is defined as 

Vl2"»=[l-<Vl3^»n-13>] Vl3 

Cn_i 2 ( 0 ) represents that part of Xn_i 3 that is eligible for promotion. 
Applying promotion and reenlistment rates to Cjj_i2(0)« part 

that is promoted and reenlisted is defined as Djj(O). 

°n<'» = 'n-l3 1i-13 Vl2<« 

Dn(0) is a part of P^{1) if this move is possible. 

Figure II-5B 

If promotion is possible from cell X^^.-y of the original in- 
ventory, then the inventory in cell Xj^.y must be subjected to 
attrition and retirement rates. That part of Xjj _7 that survives 
these rates is defined as Cjj_^(0). 

VeW’ = [l-<V7^V7> ] \-7 

Cn_^(0) represents that part of X^.y that is eligible for promotion. 
Applying promotion and reenlistment rates to Cn-6(0)> that part 



that is promoted and reenlisted is defined as Dn(0)- 
■’„«» = *n-7 ^n-7 

Dji(O) is a part of Pjj(l) if this move is possible. 
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Figure II-5C 

If the passage of one year can move the inventory in cell 
into Pn(l)« then ^ must be subjected to attrition and 
retirement rates. That part of which survives these rates 

is defined as Cjj(O). 



C^(0) represents that part of that is eligible for promotion. 

Applying promotion rates and reenlistment rates to Cj^(O), that 
part that is reenlisted but not promoted is a part of Pn(l) this 
move is possible. 

Figure II-5D 

All cells, Pj^(l)i eligible for a training input, T^(l). 

By using Figure II-4 to determine which subset of the 
possible inputs shown in Figure II-5 is applicable to each cell 

2, . . . , 72), the following set of projection equations 

are derived: 
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p„(l) =< 



D (0) + T (1) 
n n 

T (1) 



0(0) = 



where 

J Vi 

0 



D (0) = 



^n-7 ^-7 
n-13 Xi-13 



n ?^1,7,19,31,49 
n = 49 

n = 1.7,19,31 

n ^1,7,19,31,49 
n = 1,7,19,31,49 



n = 8 - 13 

n = 20-30,32-43,49,61 
n= 1-7,14-19,31,44-48,62-72 



With the projection equations defined for every cell of the 
projected inventory matrix, the next step is to compare these 
equations against the stated requirements, allowing the training 
inputs to be unknowns. This provides the necessary inputs to a 
linear program. The linear programming method that is utilized 
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is as defined in GASS, L4 J , In naatrix formulation this method 

is used to 

minimize cY , subject to AY ^ B 
where 

c = j^Cj J s Cost Coefficient Vector 
Y = j^Yj J = Unknown 
cY = ] = Objective Function 

A = J = Matrix of Coefficients 
B = j^B^ J = Restraint Vector 



and in this problem 

AY > B represents an array of equations 

generated from comparing the projection equations with the 
stated requirements, with the Y matrix representing the un- 
known training input. 

The necessary inputs to the linear program are therefore 
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the matrices A, B, c, and cY. It is therefore these matrices 
that must be derived from the projection equations and the re- 
quirements. 

The method of linear programming is dependent on the manner 
in which the requirements are stated. The problem can be 
approached assuming that the requirements could be stated for 
each of the 72 cells of the projected inventory matrix. The re- 
sultant matrices A, B, c, cY are discussed in Appendix II; 
however, this approach is not pursued further as this form of 
stating requirements is not anticipated to be realistic. This 
leads to a form wherein the requirements are stated in terms of 
requirements for entire pay grades. These requirements are 
defined as 

Rj(m) = The sum of the total requirements for each of the 
pay grades E-3 through E-9 for the m^^ year of 
projection. 

R 2 (m) = The sum of the total requirements for each of the 
pay grades E-4 through E-9 for the m*^ year of 
projection. 

R 3 (m) = The sum of the total requirements for each of the 
pay grades E-5 through E-9 for the m^^ year of 
projection 
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R 4 (ni) = The sum of the total requirements for each of the 
pay grades E-6 through E-9 for the m^^ year of 
projection. 

R 5 (m) = The sum of the total requirements for each of the 
pay grades E-7 through E-9 for the m^^ year of 
projection. 

It should be noted that this form of stating requirements 
demands that an input to the model must be the total requirement 
for each of the pay grades. The following definitions are there- 
fore in order: 



REQi(m) 
REQ2(m) 
REQ3(m) 
REQ^Cm) 
REQ5(m) 
and therefore 



Total requirement for pay grade E-3 for the 
m^^ year of projection. 

Total requirement for pay grade E-4 for the 
m^^ year of projection. 

Total requirement for pay grade E-5 for the 
m^^ year of projection. 

Total requirement for pay grade E-6 for the 
m^^ year of projection. 

Total requirement for pay grade E-7 through 
E-9 for the m^^ year of projection. 



39 



R^(m) 



R2(m) 



R3(m) 



R4(m) 



Rg(m) 



k=l 



REQj^(m) 



y REQ. (m) 

k^2 



y REQ. (m) 
k=3 



^ REQj^(m) 

1^4 



REOg(m) 



40 



•r*. -«:i* 

-li 



This procedure for the specification of requirements appears 
to be the most useful as it approaches the problem from an 
operational point of view. Merely stating the requirements for 
individual pay grades offers the following dilemma. 

Consider the requirements for pay grade E-5 to be 1000 men 
and for pay grade E-6 to be 500 men, and the projected inventory 
to be 800 men in the E-5 and 700 men in E-6. Such a situation 
would fail to meet the requirements for pay grade E-5; however, 
in actuality, since there exists an excess of 200 men in E-6, the 
requirements are really oversatisfied with E-6 personnel filling 
E-5 billets. By pyramiding the requirements, as defined by the 
Rjj(m) above, this dilemma is eliminated. Such a form of stating 
requirements will, however, necessitate consideration as to the 
varying costs of training personnel in different pay grades. A 
discussion of this cost analysis can be found in Appendix A, from 
which is generated the definition: 

= Cost coefficient associated with the pay grade 



(k = 1, 2, . , , 5) 



c 
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where kj — E - 3 




E - 4 



k- 



■3 



E - 5 



k 



‘4 



E - 6 



k 5 — E - 7 through E - 9 



With a specific form of requirements available, with 
associated cost coefficients, it is possible to derive the matrices 
A, B, c , cY, and hence the necessary inputs to a linear program 
which will compute the minimum training input necessary in order 
to meet the requirements. This derivation, with a listing of the 
appropriate solutions is shown in detail in Appendix B. 

With the training inputs now available, the projection 
equations will determine the projected inventory at the end of 
one year of projection. 
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C. Inputs and Outputs 



In this section the inputs and outputs of the basic model are 
delineated and defined, 

1. Inputs 

a. Initial Inventory 

The initial inventory is the basic group of personnel 
under consideration whose future composition is of interest to the 
personnel planner. This group is subdivided by pay grade, obligated 
service, and current enlistment into a matrix of 72 cells as shown 
in Figure II-3, In this thesis the term "initial inventory" can apply 
to any logical subdivision of the personnel system that can be 
described as "critical", in the same sense that any scarce but 
valuable commodity might be regarded as "critical". Thus the 
initial inventory may be composed of: 

(1) A critical rating (i. e. , the ET rating); 

(2) A critical NEC (i, e. , NEC 3322, Inertial Navigation 

Subsystem Technician, Mk II); 

(3) Any desired combination of ratings and/or NEC's. 

The criteria for allocating individuals of the subject 

personnel group into their appropriate rows and columns (i. e. , 
cells) of the inventory matrix are given below: 



43 



(1) Row Criterion 

This is a breakdown by obligated service. An individual 
in the inventory is represented by a row numbered (n+1) according 
to the relationship: 

n = (calendar year of expiration of active obligated service) 

- (present calendar year) 
where n can take on the values n=0, 1, 2, 3, 4, 5. 

Example: A person in a calendar year 1964 inventory whose 
expiration of active service (EAOS) occurs in calendar year 1967 
would be placed in row number four. That is, 
n = 1967 - 1964 = 3, so that 
row (n+1) = row (4). 

(2) Column Criterion 

This is a breakdown by pay grade and current 

enlistment. 
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COLUMN 


PAY GRADE 


CURRENT 

ENLISTMENT 


LENGTH OF SERVICE 
X(IN YEARS) 


1 


E-3* 


First 


0 < 


X 


<_ 6 


2 


E-4 


First 


o 

A 

X 


< 6 


3 


E-4 


Second 




X 


> 6 


4 


E-5 


First 


0 < 


X 


< 6 


5 


E-5 


Second 




X 


> 6 


6 


E-6 


First 


0 < 


X 


< 6 


7 


E-6 


Second 


6 < 


X 


< 12 


8 


E-6 


Third 




X 


> 12 


9 


E-7>!< 


Second 


6 < 


X 


< 12 


10 




Third 


12 < 


X 


< 18 


11 


E-7** 


Fourth 


18 < 


X 


< 24 


12 


E-7** 


Fifth 


24 < 


X 


< 30 


Notes: * 


Pertains to pay grades E-1, 2 and 3 


inclusive. 







siois Pertains to pay grades E-7, 8 and 9 inclusive. 

The initial inventory used in evaluating the computer 
models was obtained from an arbitrarily chosen combination of 
four specific NEC's, encompassing two separate critical ratings, 
associated with the Fleet Ballistic Missile program, and repre- 
senting closely related skills and training. The authors have 
chosen not to identify the specific NEC's and ratings used in 
order to preclude the necessity of giving the thesis a military 
security classification. The initial inventory was compiled from 
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data obtained from the Bureau of Naval Personnel, using the 
criteria enumerated previously. Figure II-6 is the standard 
initial inventory used for all computer runs discussed in this 
thesis. 



PAY GRADE 


E-3 


E-4 


E-5 


E-6 


E-7 


ENLISTM] 

Years 

Obligated 

Service 


ENT 


1 


1 2 


1 2 


1 2 3 


2 3.4 5 


1 


0 


5 


0 


13 


3 


1 


1 


0 


0 


0 


4 


0 


2 


6 


12 


0 


18 


10 


0 


11 


1 


1 


5 


7 


0 


3 


8 


33 


1 


16 


12 


4 


15 


9 


2 


9 


17 


0 


4 


3 


13 


1 


21 


13 


3 


13 


13 


0 


14 


18 


0 


5 


20 


20 


0 


32 


12 


5 


8 


3 


5 


11 


4 


0 


6 


0 


48 


0 


26 


2 


2 


7 


4 


4 


13 


2 


0 



Figure II- 6 

Standard Initial Inventory 
b. Requirements 

The Five Year Force Structure and Financial 
Program is a top secret document by which the Secretary of 
Defense establishes the approved force structure and financial 
levels for the various military departments for the succeeding 
five fiscal years. Emanating from this document are the force 
levels for each of the approved programs or weapons systems in 
the budget, including the total cost of military personnel for each 



46 




i 



program. Thus each approved program includes the total cost of 
the personnel sybsystem associated with each particular weapons 
system. The overall personnel requirements for each approved 
program are formulated within the budgetary constraints of this 
document. These requirements are then broken down within each 
service to determine the numbers of personnel of various skills, 
qualifications, and experience required. 

It is assumed in all the models described in this thesis that 
the specific requirements are known by pay grade per year for 
personnel corresponding to the same grouping as the initial in- 
ventory, It is noted that these input requirements by pay grade 
per year are different from the kind of requirements used in the 
restraint equations of the mathematical model in Chapter II, 
page 38, The pyramided form of requirements found in the re- 
straint equations are computed from the input requirements 
discussed in this section. It is these input requirements, of 
course, that are of primary interest to the personnel planner. 

Attempting to determine the requirements corresponding to 
the group of critical NEC's used for the standard initial inventory 
in Figure II-6, the authors found that the data furnished by the 
Bureau of Naval Personnel showed that these particular require- 
ments were broken down into two categories, "supervisors" and 
"technicians", where supervisors are defined to be personnel of 
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pay grades E-5 through E-9, inclusive, and technicians to be men 
in pay grades E-1 through E-4, inclusive. Upon comparing the 
relative percentage breakdown of the initial inventory by pay grade 
with the proportional breakdown of actual requirements by super- 
visors and technicians, an arbitrary percentage distribution of 
requirements by pay grade was determined. The table in Figure 
II-7 should help the reader to understand fhe derivation of this 
distribution of requirements. Note that each percentage in the 
table represents the appropriate ratio of personnel in that pay 
grade to the total of personnel in all pay grades. 



Pay 

Grade 


Initial 

Inventory 


Actual BuPers 
Requirements 


Arbitrary Require- 
ments Breakdown 


E-7 


ZO. 5% 


Supervisors 

70% 


15% 


E-6 


17. 4% 


20% 


E-5 


30. 9% 


35% 


E-4 


Z4. 8% 


Technicians 


25% 


E-3 


6.4% 


30% 


5% 



Figure II-7 

Derivation of Requirements 
Since the total requirements per year for five 



years force structure are known, the requirements by pay grade 
per year are obtained by applying the arbitrary distribution 
shown in the right hand column of Figure II-7 to the total require- 
ments per year. Figure II-8 is a tabulation of the standard 
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requirements by pay grade per year for the personnel grouping 
used in all computer runs in the thesis. 



Year 

Pay 

Grade 


1 


2 


3 


4 


5 


E-7 


76 


89 


99 


106 


106 


E-6 


102 


120 


133 


142 


142 


E-5 


177 


210 


51 


248 


248 


E-4 


128 


149 


166 


178 


178 


E-3 


26 


31 


33 


36 


36 



Figure II-8 

Standard Requirements 
c. Attrition Operator 

Associated with the passage of each increment of 
time is the likelihood that a certain fraction of the initial inventory 
will have been lost or "attrited" due to unforseen circumstances 
which may be called "acts of God. " In order to account for this 
"unplanned" loss of personnel the attrition operator ( ) is 

introduced into the model. This operator will be composed of the 
summation of specific types of losses and expressed as a percentage 
reduction of the total inventory per year. 

The numbers in parentheses following the various types 
of losses are the specific item numbers for attrition of active 
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duty enlisted personnel as defined in NAVPERS 15658, Navy and 
Marine Corps Military Personnel Statistics [ 5 ] . 



By definition the attrition operator consists of 
the following types of personnel losses: 

(1) Death (f!U:Z] 

(2) Physical Disability (#1120) 

(3) Discharge by Undesirable Discharge (#1000) 

(4) Discharge by Bad Conduct Discharge (#1040) 

(5) Discharge by Dishonorable Discharge (#1050) 

(6) Discharge by Reason of Dependency or Hard- 
ship (#890) 

(7) Intra-Navy Transfer (#1250) 

(8) Miscellaneous 

The actuarial data used to determine the attrition 
operator matrix should correspond to the same group of personnel 
in the initial inventory. In reality, it is ^recognized that this in- 
volves a degree of data specificity that does not exist. Thus the 
user of this model, as a substitute for a complete breakdown of 
actuarial data by inventory groups, must make certain value 
judgements or estimates based on the data available in order to 
apply it to the particular inventory under consideration. 

Conceptually it is possible to determine a unique 
value for each cell of the attrition operator matrix. Because of 
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the lack of discrimination of basic attrition data with regard to 
the specific sample in the initial inventory as well as the intuitive 
idea that an "act of God" may be equally likely to occur at any 
particular cell of the inventory matrix, a uniform value for the 
attrition operator was used throughout. The standard value used 
in all computer runs is O' n “ determined from 

data found in reference [ 5 j as follows: 



Category Percentage 

Death 0. 16 

Physical Disability 0. 28 

Undesirable Discharge 0. 43 

Bad Conduct Discharge 0. 19 

Dishonorable Discharge 0. 00 

Dependency or Hardship 0. 09 

Intra-Navy Discharge 0. 26 

Miscellaneous 0. 09 

Total 1. 50% 

d. Retij ement Operator 



A second means by which personnel leave the 

service is by voluntary retirement, and thus this factor must 

also be represented in the model. The retirement operator 

( P ) is defined to be the percentage of personnel represented 
n 

in the n^^ cell of the inventory matrix who are separated through 
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voluntary retirement by reason of having completed twenty years 
or more of satisfactory active service. 

The values assigned to each cell of the retirement 
operator matrix are inputs to the model and hence are at the dis- 
cretion of the user. The retirement operator used throughout 
this thesis has positive values for only eleven cells of the basic 
inventory matrix, (i. e. , cells numbered 62 through 72), because 
only these cells represent a total length of service of twenty years 
or more. The remaining 61 cells, (i. e. , numbered 1 through 61), 
all have a value of zero to represent the fact that voluntary re- 
tirement from one of these cells is not possible. 

As for the attrition operator, the statistical 
data available in the Bureau of Naval Personnel for the retirement 
operator is not sufficiently broken down by NEC or rating to be 
used directly. It is necessary, therefore, for the personnel 
planner to make an estimate of the values of the retirement operator 
matrix for use in this model. The positive values used as the 
standard for computer runs in this thesis were estimated by the 
authors, and are tabulated below. 
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Approximate 
Length of Service 
(years) 



62 


0 . 50 


20 


63 


0 . 25 


21 


64 


0 . 50 


22 


65 


0 . 50 


23 


66 


0 . 50 


24 


67 


0 . 50 


25 


68 


0 . 50 


26 


69 


0 . 50 


27 


70 


0 . 75 


28 


71 


0 . 75 


29 


72 


0 . 75 


30 



e. Promotion Operator 

In any simulation of the enlisted personnel 
system, the function of promotion, or the process whereby 
personnel are advanced in seniority level within the structure, 
must be considered. In the basic model of this thesis, the 
promotion operator n)is defined to be the ratio of the number 
of men promoted to the total number of men represented in the 
n^^ cell of the inventory matrix and is expressed as a percentage. 
The terms "promoted" or "promotion" used herein are specifical- 
ly limited to advancement from one enlisted pay grade to the next 
higher pay grade. It is noted that promotion from enlisted to 
commissioned officer status is, by definition, not considered to 
be part of the promotion operator, but is considered elsewhere in 
the model. Specifically, promotion to commissioned status is 
included in the attrition operator under the subcategory of "intra- 
Navy transfers". 
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Because the promotion operator is limited to 
intra-enlisted promotions, and the inventory matrix compresses 
pay grades E-7, E-8, and E-9 into one group, (identifying them all 
as pay grade E-7), the highest promotion considered in the model 
is that from E-6 to E-7. Thus the cells of the promotion operator 
matrix which correspond to pay grade E-7 in the inventory matrix, 
(cells 49 through 72) have a zero value, to indicate that it is not 
possible to be promoted from one of these cells. Similarly, 
because the inventory matrix compresses pay grades E-1, E-2 
andE-3 into a single group, (identifying them all as pay grade 
E-3), the lowest promotion considered by the model is that from 
E-3 to E-4. 

As was the case with the other operators, the 
promotion operator matrix should be determined from statistical 
data that corresponds to the same sample population as the initial 
inventory. While it is felt that more definitive data exists in the 
Bureau of Naval Personnel for the promotion operator than for 
any other operator, the authors estimated the values used for 
the standard promotion operator in order to expedite the develop- 
ment of the computer models. 

Figure II-9 is the standard promotion operator 
matrix used for all computer runs to evaluate the models. 
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A further distinction exists between the promotion 
operator and the other operators. From the planner's standpoint 
promotion can be thought of as a controllable" variable* whereas 
attrition, reenlistment and retirement might be thought of as 
"uncontrollable" variables. A controllable variable is one which 
the Bureau of Naval Personnel can vary by means of policy changes; 
an uncontrollable variable is one whose value is determined 
primarily by actuarial data, and over which the planner can 
exercise little or no control. 

In a limited sense, retirement and reenlistment 
might also be considered controllable variables in the personnel 
system. An example of such a situation might be a period of 
national emergency in which retirements and enlistments are 
"frozen", leaving only attrition or act of God loss as the only 
method by which a person could leave active service, 
f. Reenlistment Operator 

As each enlisted man approaches the end of 
his current enlistment he must decide whether to reenlist and 
thus continue his period of service, or to accept a discharge and 
leave the armed forces. Any model which simulates a personnel 
system should provide for this decision process. In this thesis, 
the reenlistment operator (5 n) defined to be the percentage 
of men represented in the n*"^ cell of the inventory matrix who 
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remain on active duty for an additional year of service. If the 
standard values of the reenlistment operator matrix (see 
Figure II-IO) are viewed in conjunction with the structure and 
function of the model, it can be seen that the reenlistment 
operator represents a more conventional reenlistment rate, i. e. , 
the ratio of the number of men reenlisting to those men eligible, 
upon completion of their obligated service. 
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Note that all cells in rows two through six 
contain the value of 1. 00 to indicate that 100% of the personnel in 
those cells will remain on active duty for an additional year. Only 
in row one, which represents personnel with one year of obligated 
service and who must soon face the reenlistment decision, do 
values of the promotion operator between 0 and 1 (0 <_ ^ 1 ) 

occur. The percentage of personnel in row one who reenlist and 
remain in the inventory matrix thus "move" with the passage of 
a year to a new cell in row six and in a different column repre- 
senting the appropriate next higher enlistment. In this fashion, 
then, the reenlistment operators for the cells of row one can be 
considered to be the effective reenlistment rates for the inventory. 

It is conceptually possible to have a unique 
value for each cell of the reenlistment operator matrix. This 
seems consistent with the notion that various lengths of enlist- 
ments are possible in the enlisted personnel system, such as 
six year, four year, and minority enlistments of less than four 
years. Concerning reenlistments, the authors found that for the 
whole Navy in FY 63, 63. 2% of all reenlistments were for a six 
year period while the remaining 36. 8% of the reenlistments were 
for periods of four years or less. (Only reenlistments made 
within 24 hours of discharge were considered in these figures. ) 

In order to obviate the complexities arising from trying to 
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accommodate different lengths of enlistments in the model, it 
was decided to introduce a simplifying assumption that only six 
year reenlistments would be permitted. This assumption is con- 
sidered reasonable and justifiable when it is remembered that 
most critical personnel subsystems require such lengthy periods 
of training that a man is required to have a six year service 
obligation prior to entering the program. 

As was the case with the other operators, 
the data used to determine the reenlistment operator matrix 
should correspond to the same group of personnel represented 
in the inventory matrix. The values for the standard reenlist- 
ment operator matrix were derived from data found in [ 6 J , 

2. Outputs 

a. Projected Inventory 

The projected inventory is the matrix which 
represents the group of personnel under consideration after the 
passage of one year. Basically the projected inventory consists 
of the initial inventory brought forward in time under the 
operation of the four operator matrices, together with training 
inputs added to the system in sufficient numbers to ensure that 
projected inventory meets requirements. 

The projected inventory matrix, corresponding 
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to the inventory matrix in structure, has 72 cells which represent 
personnel with specific obligated service, pay grade, and current 
enlistment. 

b. Training Input 

The training input is the total quantity of person- 
nel added to the particular personnel system under consideration 
to ensure that the projected inventory meets requirements. The 
training input is a matrix of 72 cells with the same structure as 
the initial inventory matrix. 

Within the conceptual framework of the model 
and the definition of the training input, it is possible to add men 
to the inventory through any of the 72 cells of the training input 
matrix. It was deemed feasible, however, from the standpoint 
of practicality to limit the training input to certain cells according 
to an arbitrarily assumed distribution. For further discussion 
of the training input distribution see Section 5 of Appendix A. 
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CHAPTER III 



ADVANCED MODELS 
A. Two Year Projection 

1. Concept 

With the construction of the one year model completed, 
the next step is the construction of a similar model that will 
perform the same functions over a two year period. It is obvious 
that "an answer" could be generated by merely cycling the one 
year model twice, using the solution to the first year of pro- 
jection as the input to the second cycle. The use of a two year 
projection model permits decisions as to quantities of the train- 
ing inputs during both years simultaneously. It can be shown 
that there exist situations where the requirements can be met 
at the end of the second year more economically by providing 
the necessary input in the first year as a lower pay grade, and 
hence at a lower training cost. This refinement becomes more 
important with a model which projects over an increasing number 
of years. 

The initial step in construction of the two year model is 
to develop a set of equations that will project the original inventory 
for a two year period. 

2. Formulation 

The projection equations must represent the inventory 
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of each cell of the inventory matrix at the end of a two year 
period in terms of particular quantities. These quantities are 
(a) original inventory, (b) operators, and (c) training inputs for 
each year. 

Using the projection equation as found on page 36 , and 

applying it to a two year period, the equation states that the inventory 
at the end of the second year is a function of the inventory at the 
end of the first year plus the training input during the second year. 

P(2) = g[ P(l) J + T(2) 

Similarly, the inventory at the end of the first year is a function 
of the original inventory plus the training input during the first 
year. 

P(l) = f [X] + T(l) 
therefore, 

P(2) = g{f[ X] + T(l)} + T(2) 

In this form P(2) is expressed in terms of known inputs (original 
inventory and operators) plus the unknown training inputs during 
both years. 

A complete derivation of these projection equations can be 
found in Section la of Appendix B. The basic form of the two year 
projection equations is: 
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P„(2) =< 



Vl [ Vl ] Vl<'> ^ Vl<l> ^ ^ 1,7,19.31,49 

D (1) + T (2) n = 49 

n n 

T(2) n = 1,7,19,31 



C (1) =< 



[ 1- < Vl^ Vl>][ V 2 ' Vl<°> ^ Vl'°> ^ Vl(l>] 

n = 3-6,9-18,21-30,33-48,51-72 

[i-‘Vi^W][Vi«»^Vi<w " = “ 

[l-WVl' ] Vl<l> n= 2,8,20,32 



n = 1,7,19,31,49 



D (1) 



'*n-7V7[l- <"n-7^ WlV8<l-V8> V7<°>^V7"»^V7'«] 

n = 9-12 

V 7 V7t-<%-7^ V 7 )] V7'« 

n-13 Xi-13 [^"''n-13''' ^n-13 ][®n-14* ^"\i-14* '^n-13*°* 

+ D ,,,(0)+T ,.(1)1 n = 21-30,33-42,49-60 

n-13 n-13 J 

Vl3 Vl3 [l-'"n-13" Vl3>] Vl3<« " = 

0 n = 1-7,13-19,31,43-48,61-72 
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It may be of some benefit to the reader to appreciate 



visually the possible moves into and out of each cell of the in- 
ventory array during the two year period. This information is 
shown in Figure III-l 

With the projection equations completed the next step is 
to provide the necessary inputs to a linear program. The linear 
program will minimize the training input for both years while 
meeting the requirements for both years. It is necessary then 
to construct the requirement restraint equations as was done in 
the one year model. 

The reader will recall that in the one year model, 
Chapter II, the solution to the training input, although originally 
determined by pay grade totals, needed to be distributed by some 
rule throughout the inventory matrix in order to represent the 
projected inventory cell by cell. This distribution now becomes 
of paramount importance since the projection equations operate 
on individual cells and therefore require the inventory at the end 
of each year to be expressable cell by cell. 

With an understanding that such a distribution rule is a 
necessary input to the model, the generation of the restraint 
equations and hence the necessary inputs to a linear program is 
possible. A detailed derivation and a listing of these equations 
and associated linear program inputs can be found in Section la 
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CELL NUMBER AFTER TWO YEARS 
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FIGURE ni-l TWO YEAR 

MOVEMENT GRAPH 
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of Appendix B. 



This linear program will compute the minimum training 
inputs for both years of projection that are necessary in order to 
meet the stated requirements for both years. 

With the training inputs now available, the projection 
equations will determine the projected inventory at the end of 
two years. 
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B. Five Year Projection Model 
1, Concept 

With the one and two year models completed, the frame- 
work from which a more complex model may be generated is 
available. In order to accomplish this, one must fully appreciate 
the iterative nature of the projection equations and the restraint 
equations as the number of years of projection is increased. In 
this section it will be shown how the projection equations are 
formed for the pay grade E-6. The equations for all other pay 
grades are derived similarly and hence will not be illustrated. 
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For a two year projection: 



P(2) = 



The number of personnel represented in approp- 
riate cells of the initial inventory projected for 
two years under the influence of the attrition, 
retirement and promotion operators, with no 
training input. 

That portion of the first year training input to 
E-5 that is promoted to E-6, and not lost from 
the inventory due to the retirement, attrition, 
or reenlistment operators during the second year. 
That portion of the first year training input to E-6 
that is neither promoted nor lost from the in- 
ventory due to the retirement, attrition, or 
reenlistment operators during the second year, 
e training input to E-6 during the second yea 
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Rewriting this equation in symbols, 




where, 

Yn(m) = 

k(m) = 



di 



= " ■'l " <*2 ^4'« " ^4<^> 

n= 31 

Inventory in cell n after m years, no training input 
Total training input to the pay grade during the 
year; k = 1, 2, . . . , 5 

Constant less than unity; dependent on known pro- 
motion, retirement, attrition, and reenlistment 
operators. 
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For a three year projection: 





The number of personnel represented in approp- 
riate cells of the initial inventory projected for 


Pn( 2 ) = 


three years under the influence of the attrition, 
promotion, retirement and reenlistment operators, 
with no training input. 


+ 

+ 


Those portions of d^T 3 (l), d 2 t ^(1) and T 4 (Z) 
that are neither promoted nor lost from the in- 
ventory due to the retirement, reenlistment or 
attrition operators during the third year. 

That portion of ''’ 3 ( 1 ) that was not promoted 
during the second year, but was promoted to 
E -6 during the third year and not lost from the 
inventory due to the retirement, reenlistment or 
attrition operators during either the second or 
third year. 

That portion of the second year training input to 


+ 


E-5 that is promoted to E -6 and not lost from the 
inventory due to the retirement, reenlistment or 
attrition operators during the third year. 

That portion of the first year training input to E-4 


+ 


that was promoted to E-5 in the second year, pro- 
moted to E -6 in the third year, and not lost from 
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the inventory due to the retirement, reenlistment 



or attrition operators during either the first, 
second, or third year. 



+ I The training input to E-6 during the third year, 
and rewriting this equation in symbols, 



18 48 

) P^(2) = ) Y^l 

n^l n4l 



(3) + d3d^T3(l) + d^d^T^Cl) + d3T^(2) 



dgd-di) T3(1) +d^T3(2) + dg^(l) + t^(3) 



These equations illustrate clearly the iterative nature of 
the projection equations. Since the restraint equations are derived 
directly from the projection equations, this identical iteration is 
also found in the restraint equations. The coefficients of the 
unknowns in the restraint equations are the elements of the matrix 
of coefficients, which is an input of the linear program model used 
to minimize the training input needed to meet requirements. If one 
were to write the five restraint equations for each year of projection, 
one under the other, with the like training inputs aligned under each 
other, the coefficients can be said to define the elements of a 
matrix (q x q) where q equals five times the number of years of 
projection. This matrix is the matrix of coefficients, and its form 
depends on the number of years of projection for which the model 
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is written. Therefore, there exists a different matrix of coefficients, 



A, for each of the three models discussed in this paper. 
Recall that for the one year projection model; 



A = 



1 1 
0 1 
0 0 
0 0 
0 0 



1 1 
1 1 
1 1 
1 1 
0 1 



where the five rows represent the five restraint equations and the 
five columns represent the five training inputs during the one year 
of projection. 
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Recall that for the two year projection model: 



First Year Second Year 



1 


1 


1 


1 


1 


0 


0 


0 


0 


0 




0 


1 


1 


1 


1 


0 


0 


0 


0 


0 


First 


0 


0 


1 


1 


1 


0 


0 


0 


0 


0 


> Year 


0 


0 


0 


1 


1 


0 


0 


0 


0 


0 




0 


0 


0 


0 


1 


0 


0 


0 


0 


0 




^7, 1 


^8,2 


^9, 3 


^10,4 


^11, 5 


1 


1 


1 


1 


1 




^8, 1 


^8, 2 


^9, 3 


^10,4 


^11, 5 


0 


1 


1 


1 


1 


Second 


0 


^9, 2 


^9, 3 


^10, 4 


^11, 5 


0 


0 


1 


1 


1 


> Year 


0 


0 


^10, 3 


^10, 4 


^11, 5 


0 


0 


0 


1 


1 




0 


0 


0 


^10, 5 


^11, 5 


0 


0 


0 


0 


1 




re the 


rows 


represent the 


five restraint equations for each of 





two years, and the columns represent the five training inputs for 
each of the two years. 

The important thing to notice is that although the size of 
the matrix has quadrupled, the number of elements that have to be 
computed anew represent less than one quarter of the new matrix; 
in particular, these elements are confined to the lower left five-by- 
five submatrix. This fact is true for every increase in years of 
projection; therefore, the problem of determining explicit 
equations for the elements of the matrix reduces to finding 
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equations for those elements in the first five columns no matter 
how many years the model will encompass in its projection. In 
addition, the reader will note that certain elements are repeated in 
the submatrix because of the iterative nature of the projection and 
restraint equations. 

Since the complexity of the problem increases with an 
increase in the number of years of projection, a judicial choice must 
be made as to the optimal number of years for which the model 
should be designed. The usefulness of the model is limited to the 
number of years for which the requirements are stated. The Force 
Requirements are currently precisely stated for five years and 
hence the optimal model is a five year model. 

2. Formulation 

The complete derivation of the necessary inputs to the 
linear program for the five year model is of such length and detail 
that it is not shown in this paper. The explicit expressions for each 
input in terms of the known operators and original inventory can be 
found in Section 2b of Appendix B. These expressions were derived 
using the same training input distribution as was used in the basic 
model. 

The output of the linear program consists of the minimum 
training inputs over the five years that will meet requirements for 
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each of the five years. With these training inputs available, the 
projection equations can be used to express the projected inventory 
at the end of each of the five years. 
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C. Inputs - Outputs 
1. Inputs 

The following inputs of the advanced models are defined 
the same as in the basic model (see Section C, Chapter II) and 
hence will not be discussed in this section: 

a. Initial Inventory 

b. Requirements 

c. Attrition Operator 

d. Retirement Operator 

e. Promition Operator 

f. Reenlistment Operator 

There are two additional inputs to the advanced models 
which were not defined in the basic model. The first of these 
two inputs is the training input distribution, which by definition 
is the arbitrary designation of those cells in the inventory matrix 
which can receive a training input. Thus the training input 
distribution determines which cells of the training input matrix will 
have positive values and which ones will have a zero value. The 
reader will note that the linear programming solution determines 
the values of the training inputs while the training input distri- 
bution determines the location of these values in the inventory 
matrix. 
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Implicit in the definition of training input distribution 



is the assumption that the total training input by pay grade per 
year ( t |j(m)) will be divided uniformly among the designated 
cells of each particular pay grade. 

The training input distribution used for the computer 
runs of the advanced models are discussed in Section 5 , 

Appendix A. 

The second of the additional inputs to the advanced models 
are the cost coefficients. As defined in this thesis, the cost 
coefficients are algebraic factors which represent costs associa- 
ted with the training of personnel of different pay grades. These 
cost coefficients appear in the statement of the objective function 
of the linear program, which determines the optimal training 
input with regard to minimum total cost of training. The objective 
function for the advanced models is to 



minimize 




5 



where j = 2 for the two year model 



j = 5 for the five year model 

C]^ = cost coefficient for the pay grade 

T ]^(rn) = total training input to the pay grade in the 



m^^ year 
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The cost coefficients provide the planner with a means 
for quantitatively expressing the complex relationships among 
training costs for personnel of various pay grades. The cost 
coefficients may take on any appropriate units, and thus are not 
restricted to the conventional units of "dollars" normally associa- 
ted with the traditional concept of cost. For example, the planner 
might desire to let the cost coefficients represent a measure of 
the difficulty of procuring personnel for a specific training 
program, with appropriate units. 

For the computer runs of the advanced models presented 
in the thesis, it was desired to obtain cost coefficients which were 
relatively simple in concept and yet quantifiable in some sense in 
order that some significance could be attached to the linear pro- 
gramming solution. Under the assumptions that (1) the fixed over- 
head costs of operating a Navy school are constant regardless of 
who the students are, and (2) that personnel in training are in a 
so-called non-productive status relative to the rest of the Navy, 
a possible measure of the relative training costs per pay grade 
would be the ratios of the various trainees' pay scales. For ex- 
ample, it costs more to train a chief petty officer than it does a 
seaman in a given unit of time, because the chief's basic pay is 
grater than the seaman's while both are in the "non-productive" 
training status. 



79 



The standard cost coefficients used in the computer runs 
of the advanced models were derived by taking the ratio of the basic 
pay of each enlisted pay grade to the basic pay of a man in pay 
grade E-3. The basic pay corresponding to the mean length of 
service for each pay grade was the criterion used to select specific 
entries from the enlisted pay table. The cost coefficients per pay 
grade determined in this manner are listed below: 

Pay Grade E-3 c^ = 1. 00 

Pay Grade E-4 c^ = 1. 32 

Pay Grade E-5 C 3 = 1. 58 

Pay Grade E -6 C 4 = 1. 87 

Pay Grade E-7 C 5 = 2. 32 

2. Outputs 

The following outputs of the advanced models are defined 
the same as in the basic model (see Section C, Chapter II) and are: 

a. Projected Inventory 

b. Training Input 
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CHAPTER IV 



DIGITAL COMPUTER MODELS 

A. General 

The digital computer programs of the mathematical models 
presented in Chapters II and III were written in FORTRAN 60 com- 
puter language, and all computer runs were made on the U. S. Naval 
Postgraduate School's Control Data Corporation (CDC) 1604 high-speed 
digital computer. All of the programming for these models was done 
by the authors except for the linear programming subroutines (i. e. , 
subroutines "LINEAR" and "SIMPLEX") which are standard routines 
on the library tape of the Naval Postgraduate School Computer 
Center . 

The manner in which the computer programs were written was 
such as to insure that the following objectives are satisfied: 

1. To provide programs that are user oriented. The pro- 
grams are designed to provide maximum flexibility and adaptability 
for the user, and to facilitate his understanding and control of the 
input data as well as his analysis of the resultant outputs. 

2. To provide within each model the capability of projecting 
an initial inventory and computing the least expensive training inputs 
for a period of five years. 

3. To demonstrate the capabilities of the computer models 
through the use of actual data inputs. 
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4. To write the programs using notation that is consistent, 
insofar as practicable, with the notation and symbology used in the 
development of the mathematical models, 

5. To write the programs using subroutines, or separate 
"packages" which accomplish certain specific operations or functions 
common to all models. 

The authors are satisfied that each of these objectives has been 
fulfilled. 

To distinguish the computer programs for the convenience of 
the reader, the following designations will apply to the models 
previously discussed: 

Model Designation 



ONE-YEAR PROJECTION 
TWO-YEAR PROJECTION 
FIVE-YEAR PROJECTION 



PIPE ONE 
PIPE TWO 
PIPE FIVE 



Program listings for PIPE ONE and PIPE FIVE are found in 
Appendix E. The listing for PIPE TWO is not included in this thesis 
because PIPE TWO was written only to facilitate the development of 
the five-year model, PIPE FIVE. If only a two-year projection run 
is desired, of course, the user need only use the first and second 
year of the five-year solution as generated by the program for 
PIPE FIVE. 

In the remaining sections of this chapter, functional flow 
diagrams of the individual programs will be presented in order to 
provide the reader with the underlying logic used in the development 
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of the respective computer programs. The computer runs used to 
illustrate these programs were made with the standard input values 
(discussed in Section C, Chapter I and Section C, Chapter II) in order 
to provide a common basis for a quantitative analysis of the solutions 
as generated by PIPE ONE and PIPE FIVE. 

B. FUNCTIONAL FLOW DIAGRAMS 

PIPE ONE and PIPE FIVE are constructed from a series of 
subroutines. The function of these individual subroutines will be 
discussed in Section C of this chapter. The actual program listing 
of these programs may be found in APPENDIX E. 

The following illustrations are functional flow diagrams of 
PIPE ONE and PIPE FIVE: 
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m 




PIPE ONE 



4600 



READ INPUT 
DATA 




$ 




P(M, N) = CK(M, N) + T(M, 






3ALL TOTAL 



M, P,SP 



/CALL TOTAL > 
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PIPE FIVE 
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1 -i 



c. 



SUBROUTINES 



Each of the subroutines used in building PIPE ONE and PIPE 
FIVE are discussed in the following paragraphs. 

1. HOUSE ONE: This subroutine is used only with PIPE ONE. 

Its function is to compute those constants and other values which are 
needed in the iterative routines of the program but do not themselves 
change during a projection. Specifically this routine: 

a. automatically sets a zero value for BETA, re- 
tirement operator, in all cells below 62; 

b. sets a value of one in all cells of DELTA, re- 
enlistment operator, except for the cells in row 
one; 

c. converts requirements by pay grade into the 
comulative requirements form used in the restraint 
equations; 

d. calls the subroutine FACTOR. 

2. HOUSE FIVE: This subroutine is used only with PIPE 

FIVE and performs the same functions as HOUSE ONE and in 
addition: 

a. computes the constants EAl, UA, VA, and WA; 

b. sets the initial value of each cell of the A-matrix 
to zero; 
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c. sets the values of the cost coefficients into the 
first row of the A-matrix where they form the 
objective function to be used in subroutine LINEAR; 

d. calls the subroutine AMAT. 

3. AMAT: This subroutine produces the values of the 

elements of the A-matrix of the linear program. It makes use of 
subroutines SUM and FIX to solve for the values of the ay found in 
the mathematical formulation Appendix B, section 2, "Advance 
Models. " 

4. FACTOR: This subroutine computes for each pay grade 

the divisor for the training input so that the computed training input 
may be distributed uniformly over the cells designated for training 
input. 

5. VENT: This subroutine solves the basic one year pro- 

jection equations found in Chapter II. It is the basic projection 
subroutine of both PIPE ONE and PIPE FIVE. In both programs, 
VENT is used in an iterative manner and is cycled as many times 
as there are years of projection. It uses as a starting inventory the 
end-of-year inventory of the year before. 

6. BVECTOR: This subroutine produces the B-vector for 

the linear program. It solves the following equation for Bi: 

Bi = I^m k ■ (that portion of the initial inventory still on board 
in the mth year corresponding to the kth pay grade) 
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7. LINEAR and SIMPLEX: These routines solve the linear 

programming problem for the optimum training input by pay grade 
and compute the cost of this input. 

8. TRAIN ONE: This subroutine computes the training 

input for PIPE ONE and then distributes this input to the appropriate 
cells of the inventory. 

9. TRAIN FIVE: This subroutine takes the training input 

by pay grade as computed by LINEAR and distributes it according 
to the training input distribution to the appropriate cells of the 
training input matrix. (VENT is then called once again and the 
training input is added in the proper cell at the appropriate time in 
the projection. ) 

10. TOTAL: This subroutine provides a means of adding 

the cells of the inventory matrix and training matrix for use in the 
subroutine OUTPUT. 

11. EXCESS: This subroutine computes any excesses and 

shortages of requirements for use by subroutine OUTPUT. 

12. OUTPUT: This subroutine causes the output to be 

printed in the form shown in Appendix D. 
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APPENDIX A 



ASSUMPTIONS 

This Appendix discusses eight major assumptions made by 
the authors in the development of PIPE. 



applicable are critical subsystems of the enlisted per sonnel 
structure of the Navy. This assumption was made because the 
critical subsystems presently require the immediate consideration 
of the personnel planner, and these subsystems permit a simpler 
formulation of the inventory matrix. 

2. The selection of which enlistments are to be represented 
in each pay grade of the inventory matrix is based upon the assumption 
that the inventories considered be those of critical personnel sub- 
systems of the Navy. An analysis of the current inventories of a 
number of these subsystems indicated that the numbers of personnel 
in the other than the following categories were negligible: 



1. It is assumed that the inventories to which PIPE is 



Pay Grade 



Enlistment 



E-3 



1st 



E-4 



1st, 2nd 



E-5 



1st, 2nd 



E-6 



1st, 2nd, 3rd 



E-7 



2nd, 3rd, 4th, 5th 
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The inventory matrix was therefore restricted to represent 
only those personnel in the listed categories. 

3. The personnel requirements by enlisted pay grade per 
year of projection are assumed to be available to the user. Although 
the requirements in this form are inputs to PIPE, the requirements 
used in the restraint equations of the linear program are various 
sums of these input requirements. This is done so as to insure 
that the requirements for the senior pay grades are filled before 
those of the junior pay grades. 

4. Since PIPE projects critical personnel subsystems, it is 
assumed that there is an unlimited supply of personnel to fill the 
requirements of these subsystems. This means that no matter how 
large a training input PIPE determines as necessary, the priority 
enjoyed by the subsystem assures that these personnel will be 
provided from the balance of the personnel system. This is 
admittedly a program of suboptimization within the personnel 
system of the Navy; however, it does appear to be a realistic 
approximation of current policy. 

5. The construction of the inventory matrix requires the 
assumption of a method by which the training input will be distribu- 
ted among the cells represented by any one pay grade. PIPE has 
been constructed to distribute the training input uniformly among 
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a pre-selected number of cells within each pay grade. The particular 



distribution selected by the authors is illustrated in Figure A-1. 
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□ = no input 



FIGURE A-1 

STANDARD TRAINING INPUT DISTRIBUTION 



The reasons for the omission of some cells are: 

a. Cell numbers (1, 7, 13, 19, 25, 31, 37, 43, 49, 55, 61, 67). 
These cells represent the first year of an enlistment and this would 

be the year spent in schools preparing the personnel for input into 
the program. 

b. Cell numbers (6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 72). 
These cells represent the last year of an enlistment and it is doubtful 
that personnel would be sent to school so late in an enlistment that 

they had but one year obligated service upon return to the operating 
forces. 
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c. Cell numbers (8, 17, 32-35, 62-65, 68-71). 

These cells are considered very unlikely to have training. 

PIPE is designed so as to provide the user with the ability to 
select the cells into which training input is to be permitted. The 
only restriction to the user's choice of distribution is that PIPE will 
automatically distribute the training input to each pay grade uniformly 
among the cells selected by the user. 

6. The manner in which training inputs are used to con- 

struct the objective function of the linear program demands con- 
sideration of the fact that the cost to the Navy is different for the 
procurement of personnel from the senior pay grades than from the 
lower pay grades. This implies that some type of cost coefficient 
must be used in the objective function of the linear program. It is 
also recognized that although this difference in cost is important 
and must be considered, there was insufficient time available to the 
authors to analyze any data that may be recorded on this subject. 
Therefore, to include this difference in cost, PIPE uses cost 
coefficients based solely on relative enlisted pay scales. The 
manner in which these cost coefficients are computed is as follows: 
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a. Compute the average base pay for each of the five 
pay grades of the inventory array 



average base pay for E-3 = PE3 

average base pay for E-4 = PE4 

average base pay for E-5 = PE5 

average base pay for E -6 = PE 6 

average base pay for E-7, 8 , 9 = PE79 



b. Divide each average pay by the average pay of the 

pay grade E-3. 

Cl = PE3 = 1 
PE3 

Co = PE4 
PE3 

C 3 = PE5 
PE3 

C 4 = PE 6 
PE3 

C 5 = PE79 
PE3 

where 

Cl = cost coefficients i = 1, 2, ... 5 

Although these coefficients do not fully express the difference 
in procurement costs among pay grades, they are considered to be 
satisfactory until such time as a cost analysis is conducted and more 
sophisticated data is available. 
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The cost coefficients are assumed to be the same for all years 
included in the inventory projection. Since these coefficients are 
inputs to PIPE, the user may change their values in any manner he 
chooses. 

7. The values assigned to the promotion, retirement, 
attrition, and reenlistment operators used in PIPE are assumed to 
be invariant during the period of projection. If the user desired to 
vary these values the symbols representing these operators must 
be additionally subscripted to provide for the designation by year, 
and the mathematical equations altered to reflect this sophistication. 

8. It is assumed that all personnel that reenlist do so for 
six years, and the time at which they reenlist is restricted to the 
normal expiration of current enlistment. Therefore neither early 
nor broken service reenlistments are considered in PIPE. 
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APPENDIX B 



MATHEMATICAL FORMULATION 
I BASIC MODELS 

a , Derivation of Inputs to the Linear Program 

1 . Requirements Specified for Every Cell 

Since cx, p , y , 6 , and X are all inputs to the model 
and hence known quantities for n = 1 ,2 , » . . ,72 , it follows that 

(1) P^(l) = T^(l) + CK^(l) n = l,2,,,,,72 

where CK^(l) is a function of only the known quantities 
enumerated above. 

The requirements state that 

(2) P^(l) > R^(l) n = 1,2 72 

where R^(l) = requirement for cell n for the first year of projection. 



(3) 


.'. T (1) + CK (1) > R (1) 
n n n 


n = 1,2, . , 


. . ,72 




or T (1) > R (1) - CK (1) 
n n n 


n = 1,2, , < 


, . ,72 


(4) 


Let B (1) = R (1) - CK (1) 
n n n 


n = 1,2, . . 


,0,72 


(5) 


IV 


n = 1 , 2 , 0 . 


, . ,72 



This yields a set of 72 independent equations, and hence 
a 72 by 72 A matrix, a 72 by 1 B matrix, a 1 by 72 c matrix, and a 
72 by 72 cY matrix. This procedure is not pursued further. 



98 




^1^ III 1 1 • • )» 



2o 



Requirements Specified for Each Pay Grade 



( 6 ) 



P„(l) = T (1) + CK (1) 



The requirements state that 
72 



(7a) 

(7b) 



y p„(i) > Rj(i) 

n=l 

J%n(l) > Rzd) 



(7c) 


IV 


Rad) 






(7d) 


y p (1) > 

n=31 


R4d) 
















(7e) 


y p„(i) > 

ni49^ 


Rjd) 






where the 


R^(l), (k= 1,2,. 


0 0,5) are as defined in 


Substituting equation (6) into equations ( 


7a= 


-7e) , we 




72 


72 






(8a) 


y T (1) > Rji) 

n4i ^ 1 


- y CKn(l) 
n=l 


= 






72 


72 






(8b) 


y T (1) > Rod) 

n^7 ^ 2 




= 






72 


72 






(8 c) 


y T (1) > R3(1) 

n=19 


- y CK (1) 

n^l9 ^ 




®4 




72 


72 






(8d) 


y Id) > R^d) 
n=3 1 


- y CK (1) 

n=3l 


~ 






72 


72 






(8e) 


y T ( 1 ) > Rjd) 

n=49 


- 

n=49 


= 


«6 



2,. _ ,72 



Chapter II. 
have 
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CK (1) = 



where the ( i = 1 ,2 , . . , 5 ) can be computed since 
0 n = 1,7,19,31 

Vl n = 2-6,14-18,44-48,62-72 

Vlt* Vl][l-<“n-l " Pn-l)]Vl^ ‘n-7V7[l-'"n-7^ V7>] V7 

n = 8-13 

^n-1 t“\i-l] t’^'n-l'^ ^n-l^l^-l'*' n-13 \i-l3[^'^ “n-13 

?1-I30^n-13 n = 20-30,32-43,49-61 

^n-13 ^n-13 t'^‘n-13'^ ^n-13^] \-l3 n - 49 

Define t (m) = J where =the total training input for the 

pay grade during the year of projection. 



(9a) 


T,(l) = ) T (1) 
1 n^i" 


(9b) 


c 

II 

CM 




30 


(9 c) 


II 




48 


(9d) 


T^d) = 1 T (1) 








72 


(9e) 


T_(l) = y T^(l) 

^ n^49^ 


Therefore 


72 5 


(10a) 


II 

It-- 




72 5 


(10b) 


rk? ^ ^2^ 



100 



f mftKd r • ^ 



fr«k - -4 




-tT *•» 

-•i f 



(10c) 



(lOd) 

(lOe) 



y T (1) = Yt (1) 

-I'*'” 



n=49 

Therefore the linear restraint equations are: 

(lla) T^(l) + "^(l) + Tgd) + T^(l) + Tg(l) > 



(lib) 

(lie) 

(lid) 

(lie) 



r^{l) + T^d) + T^d) + Tgd) > 



T^d) + "^(i) + Tgd) > B^ 



T^d) + TjCl) > Bj 






These five linear restraint equations produce the following 
necessary inputs to the linear program: 



1 1 



1111 



0 0 
0 0 
0 0 



1 1 



1 1 
0 1 



Bo 
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c = ( , Cg , c^, Cg ) as defined in Chapter II. 

Y = t(I) = [t^(D ] 
and the objective function is 

5 

minimize cY = inin \ 

The output of the linear program is the matrix 
T (l) , the elements of which are the total training input for each 
of the five pay grades „ These totals represent the minimum training 
input to each of thepay grades that will produce a projected inven- 
tory that satisfies the stated requirements at the end of one year of 
proj action . 

Recall, T^(l) = y T^(l) 

In order to represent each cell of the projected inventory, P^(l) » it 

is necessary to determine T (1) explicitly for each n = 1 ,2, . ,72. 
n 

This demands some rule by which the total training input to any one 

pay grade, distributed among the cells within that pay 

grade. A discussion of this problem and a suitable solution is explained 

in detail in Appendix A . Using this rule, each T^(l) (n = 1,2, , . , ,72) 

can be determined once having solved for the (k = 1 ,2 , „ . „5 ) . 

Knowing the value for each T (l) , the projection equations, (the only 
n 

unknowns of which are the T (l) ) , can then be used to explicitly state 
n 

each of the P (l) and hence the entire projected inventory after one 
n 

year, P(l) » 
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APPENDIX B 



MATHEMATICAL FORMULATION 
II „ ADVANCED MODELS 
a. Two Year Model 

This section of Appendix B contains the mathematical formulation 
of the two year model discussed in Section A of Chapter III. 

1. Derivation of Projection Equations 

Applying the projection equation of the basic model 
(Chapter II, page 36) to a two year period, the inventory at the end of the 
second year of projection is expressed as; 



and T^(2) . In order to express Pj^(2) as a function of known quantities and 

the unknown training inputs, it is necessary that Pj^(2) be expressed as a 

function of C (0) , D (0) , T (1), and T (2) . The projection equations for the 
n n n n 

basic model can again be used, but now to express C (l) and D (l) in the 

n n 

desired terms as follows: 



(12) P^(2) = 




n?^l,7,19,31,49 



n=49 



n=l,7,l9,31 



In these equations P^(2) is a function of C^(l) , D^(l) 



First, C^(l) can be expressed as 




n?^l,7,l9,31,49 



0 



n=l,7,l9,3l,49 
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where 



( 14 ) 



(15) D^(l) = 



?.-2[l-\-2]VlW n^2.8,20,32,50 

D (0) + T (1) n-20"30, 32-43, 49-61 

n-1 n-1' 

T ,(1) n=l-7, 14-19, 31, 44-48, 62-72 

n-1 

Second, D^(l) can be expressed as 

^n-13Vl3 Vl2<« n=20-30,32-43, 49-61 

0 n=l-7, 14-19,31,44-48,62-72 



where we express C _(l) and C ,„(l) as 
^ n-6 n-12 



(16) C^.g(l) 



and 



(17) C„.^3(l) = 



I ['-<V 7 ^ Wh- 7 <« n^ 7 , 13 , 25 , 37.55 

I 0 n= 7 , 13 , 25 , 37,55 

j['-'Vl 3 ^®n- 13 >]Vl 3 <'> 13 . 19 , 31 , 43,61 

0 n= 13 , 19 , 31 , 43,61 



and then express P ^(1) and P ,„(l) as 
n-7 n-13 



(18) P^_^(l) = 



and 



(19) P„.j3(l) = 



'n- 8 [l-V 8 k- 7 ®> ^ ^ V 7 (l) "^ 8 , 14 , 26 , 38,5 

n=56 

n=8,l4,26,38 



V7<‘»^V7'l) 

^n-7<l> 



Vl4t-Vl4]Vl3<°! ^°„-13®> ^ Vl3'« 

n= 62 



Vl3<°> ^ Vl3<l) 
Vl3<'> 



n= 14 , 20 , 32,44 
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r 



A graphical representation of this breakdown is: 




(13) 



I 

(14) 




i I 



(18) (19) 



Equations (13) through (19) serve to express the quantities 

C (1) and D (l) of equation (12) in terms of known inputs and training inputs, 
n n 

Retracing these steps so as to ultimately be able to write equation (12) in 

the desired terms, the procedure is: 

Substitute equation (19) into equation (17) which yields 
/ 

^ Vl3>][Vl4<l-Vl4> ^ ^ Vl3<«] 

n?^ 13-14,19-20,31-32,43-44,61-62 

^ »n-13>] [°n-13™ ^ ^n-13'«] "= 

['-'"„-13"W][Vl3<«] n= 14,20,32,44 

. n= 13,19,31,43,61 



(17a) 






Then substitute equation (18) into equation (16) which yields 

^ Pn-7>][V-8(l-^-8)°„-7») ^ Dn-7® ^Tn-7<«] 
n?^ 7, 8, 13, 14, 25, 26, 37, 38, 55, 56 

[l-'V7^W]Wl> n = 8,14,26,38 

0 n = 7,13,25,37,55 



(16a) 
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Then substitute equations (l6a) and (17a) into (15) which yields 



*n-13Vl3[l-<Vl3 ^ l>-13>][Vl4<l-Vl4>Vl3® * °n-13® 
+ Tn_i 3 (l)] n = 21-30,33-42,49-60 

(15a) n-l3^n-13[^“^^n-13 ^n-13^]^n-13^^^ n = 20,32 

D„(l)= { ‘n-7\-7[l-'“n-7 ^ ^n-7> ] &n-8'l-\,-8>^n-7® ^ °n-7® ^ ^n-7<« 

n = 9-12 

V7V7[l-'“„-7^^n-7)]V7<'> "-8 

^ 0 n = 1-7,13-19,31,43-48,61-72 

Then substitute equation (14) into equation (13) which yields 

' ['-<Vl ^ ^ Vl™ ^ Vl<«] 

(13a) n 7^2,8,20,32,50 

C„(l)= I ^ " = =« 

ri-(» n+P JITt ,(1)1 n = 2,8,20,32 

L n-1 n-1 _i n-1 J 



Equations (l3a) and (15a) along with equation (l2) are therefore 
the projection equations for the two year model. Given an initial inventory 
and operator matrices, these projection equations will describe the inven- 
tory at the end of the second year of a two year projection, in terms of the 
training input during both years . 
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2 . Derivation of Inputs to the Linear Program 



Stating the requirements in the same manner as in equation 
(7) , (page 99 ) , but for both years of projection: 



(20a) 


A| 


(20b) 


> R^d) 


(20 c) 


h"' - 


(20d) 




(20e) 


n=49 


(20f) 


y p„(2) > Rj(2) 
n=l 


(20g) 




(20h) 


72 

y P„(2) > R.(2) 
nfel9" 3 


(20i) 


Pi 

A1 


(20j) 


f P„(2) > R3(2) 
n=49 
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^'1 
<* .Mof 

...<P 



The restraint equations generated from equations (20a) 



through (20 e) are identical to the five restraint equations of the one year 
model, equations (lla) through (lie) „ The problem reduces then to the 
derivation of the five additional restraint equations represented by equations 
(20Q through (20j) . In order to facilitate the derivation the following com- 
binations of operators and original inventory are defined: 




G2 = 6 
n r 



6 V 

n-13 Vi-13 









M 



n 





V 



n 




W = G2 E3 
n n n 



Y 



n 
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The training inputs are defined as in equation (9) but 



for both years of projection. 



(21a) 



(21b) 



(2lc) 



(2 Id) 



(2le) 



T (m) = ) T (m) 



n=l 

30 

"3<“> = I 

n=l9 

48 

n=3 1 
72 

T (m) = y T (m) 
5 4_i.„ n 

n=49 



m = 1 ,2 



m = 1,2 



m = 1 ,2 



m = 1,2 



m = 1 ,2 



It is at this point that some rule must be assumed for 
the distribution of the training input throughout each pay grade. The distrib- 
ution assumed is discussed in Appendix A and generates the following 
definitions; 



(22a) 


T^(m) = 1/4 T^(m) 


m = 1 ,2 


n = 2-5 


(22b) 


T (m) = l/6 T (m) 
n 2 


m = 1,2; 


n = 9-11,14-16 


(22c) 


T (m) = 1/8 T (m) 
n 3 


m = 1,2; 


n = 20-23,26-29 


(22d) 


T^(m) = 1/8 T^(m) 


m = 1,2; 


n = 38-41,44-47 


(22e) 


T (m) = 1/8 T (m) 
n 5 


m = 1,2; 


n = 50-53,56-58 


(22f) 


T (m) = 0 m = 1 , 2; n = 
n 


1,6-8,12- 


13,17-19,24-25,30 



48-49,54-55,60-72 
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With these definitions, the projection equations (12) , 



(13a) and (l5a) can be rewritten as follows: 



P„(2) 

Pn(2) 

P_^(2) 

P^(2) 

P„(2) 

Pn(2) 

P„(2) 

P^(2) 

P„(2) 

Pn(2) 

Pn(2) 

P„(2) 

Pn(2) . 

P„(2) 

P„(2) 

P„(2) 

Pn(2) 

Pn(2) 

P„(2) 

P„(2) 



1/4 t^(2) 



n = 1,7,8,19,31,32 
n = 2 

U fcKl + 1/4 T (1) 1 + 1/4 T (2) n = 3-5 

nL n 1 J 1 

Un[CKln + l/4 t^( 1)] n = 6 

U^CKl + V^rCK2 + 1/4 T (1) "I + 1/6 t.(2) n = 9 

n n nL n i j 2 

^n[^^^n n = 10,11 

UnfcKln + I/6 T2(1)] + V^[cK2^ + l/4 t^(1)] n = 12 

U CKl + 1/6 t_(2) n = 14 

n n 2 

Un[cKln + I/6 r^d)] + I/6 r^{2) n = 15,16 

UnicKln+l/ST^d)] n = 17 

n = 13,18,43,61-72 
n = 20 

U Fcki + 1/8 T-(l)1 + W CK3 + 1/8 t.(2) n = 2l 

nL n 3 _1 n n 3 

Un[CKln + l/8 T^d)] + ^^[^^3^ + l/6 T^d)] + l/8 -^(2) n = 22,23, 



U CKl 
n n 



1/8 T3(2) 



= u Fcki + i/8 t.(i)1 + w Fck 3 + i/6 t.(i)1 

nL n 3 _1 nL n « _1 



= U CKl + W CK3 
n n n n 



= U CKl + W CK3 + 1/8 tM) 
n n n n 3 

= Un[cKl^+l/8T3(l)] +W_^CK3^ 

= U CKl +wFcK3 + 1/8 t. (1)1 
n n nL n 3 J 

= U CKl + W CK3 + 1/8 tA2) 
n n n n 4' ' 



n = 24 

n = 25,37,55 
n = 26 
n = 30 
n = 33-36 
n = 38 



no 



n = 39-41 



Pj^(2) = + 1/8 t^(D] + W^[CK3^ + l/8 t^Ci)] + l/8 x^(2) 

Pn(2) = U^[CK1_^ + 1/8 t^(d] + W^[cK3^ + 1/8 Tjd)] 



P (2) = U CKl + 1/8 t.( 2) n = 44 

n n n 4 

P (2) = U [CKl + 1/8 tJD] + 1/8 t.( 2) n = 45-47 

n nL n 4 j 4 

P (2) = U [cKl + 1/8 T.(l)1 n = 48 

n nL n 4 J 

P (2) = W CK3 n = 49 

n n n 

P (2) = U D , (0) + W CK3 + 1/8 t^( 2) n = 50 

n n n-1 n n 5 



Pn(2) = Uj^[cKl^ + 1/8 Tjd)] + W_^[CK3^ + 1/8 t^(1)] + 1/8 -^(2) 

P (2) =U fcKl + 1/8 Tdl)l + W fcK3 + 1/8 t.( 1)1 
n nLn 5JnLn 4J 

P (2) = U CKl + W CK3 + 1/8 t^( 2) n = 56 

n n n n n 5 



n = 42 



n = 51-53, 
57-59 



n = 54,60 



With an equation now available for each cell of the inventory 
matrix at the end of the two years of projection, it is now possible to gen- 
erate the desired restraint equations . Because of the manner in which re- 
quirements are stated the remaining five restraint equations are more easily 
derived in reverse order. 

Recall the tenth requirement equation is 



(20j) y P^(2) > R^(2) 

n^49 ^ 

and by proper substitution from equation (23) , the tenth restraint equation is 



(24) 






■ n=51 n=57 



n=51 n=57 



72 60 

• R J2) - r y u CKl +y WCK3 +U_D..(0)1 

-5 L L-^ n n Zj._ n n 50 49 J 
n=51 n=49 
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where the appropriate elements (a^^) of the matrix of coefficients, A, and 
the elements (B^ of the restraint vector B are defined as follows; 



11,4 



n«51 n=57 



54 60 



72 60 

®n = ^5® - [ Lyn^'^ln ^ ' ®10 



n=51 



n=49 



Therefore we write equation (24) as 

^ 11,4 -^^ 1,5 V ^’ - ®11 

Recall the ninth requirement equation is 



72 48 72 

(20i) > P^(2) =) P^(2) +y P^{2) > R.(2) 



and by proper substitution from equations (23) and (24a) the ninth restraint 
equation is 

36 42 42 48 

(25) l/8[ y W„ 47 wJt3(1) 4 l/8[ 7 U 4y 4 r^(l) 

n^3 n^9 n=39 n=45 



11 



48 HZ. 

-feCl) + t^(2) + 43 ( 2 ) > R^(2) - [S 30 +7 U„CKl„4y W CK3„ 

n=33 



42 

+y ^ 

*'n=33 ' 



Defining (a^^) and (B^) similarly we have 
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42 48 

n=39 n=45 

48 42 

«10 = - [ s ,0 = '^4<2> ' [^10 ^ ^ s ] 

n^3 n=o3 



Therefore we rewrite equation (25) as 

^10,3 "3'1> -"^0,4 "4<1’ -"^1,5 Vl> - ®10 



Recall the eighth requirement equation is 
72 



(20h) 



72 

Zi 

n=l9“ 



Z,^n<2) =? Pn<2) tSl""*'' - 

n=19 



(2) 



and by proper substitution from equations (2 3) and (25a) the eighth restraint 
equation becomes 



(26) 



24 

‘"• LL '’" 



29 24 30 

ri“X / n — n — 27 



10,3 



^(1) 



i '^3<'> - [^10 ^ )] 

Defining (a_) and (B^) similarly, we have 

a- - = l/6[ y W +\ W 1 

'-n^2 " n^27 
24 30 

30 

= '^ 3 ® - [®10 ^ ®9 ^ )] = '^ 3 <'> ' [® 10 ^ ^ 9 ^ ® s ] 

n=2 1 
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Therefore we rewrite equation (26) as 






+ ^5(2) ^ Bg 



Recall that the seventh requirement equation is 

E l§ 72 

P (2) =y P (2) +Y P (2) > R (2) 

and by proper substitution from equations (23) and (26a) , the seventh restraint 
equation becomes 

(27) l/4[ f Vj^(l) + l/6[ fu U ] T^d) + a r^Cl) + Tjd) 

n=9 n=10 n=15 



4^4*1* ■"’■2*^* ■"Y® ■^''4® ■"’'5® 

n=9 n==y 



Defining (a.^) and (B^) similarly, we have 






n=lO n^5 

«8 = - [^ 10 ^ ® 9 ^ “ 

= ■'a® - [^10^ ^8^ ®7] 



114 



Therefore we rewrite equation (27) as 



(27a) 



^^8,2"2''> 

+ j^{2) + t^{2) + t^( 2) + ”^(2) ^ Bg 



Recall the sixth requirement equation is 

( 20 f) ?P„( 2 )=fp„( 2 )+IP„( 2 ) > Ri( 2 ) 

n=l n=l n=7 

and by proper substitution from equations (23) and (27a) the sixth restraint equation 
equation becomes 

(28) l/4[luj^(l) ■fa 3 _jT^(l) +^ 8 , 2 ^ 2 ''> 

n=3 



+ aii gTg(l) +t^(2) +t^{2) +Tg(2) +t^(2) + t^(2) > 

> Ri(2)-[Sj^^S3+S3^S,+ f uCKlJ 
n=3 

Defining (a^^) and (B^ similarly we have 

^^8,1 

n=3 

®7 = V ®7^Zyn°"'n] = - [®10^ ®8^ ®7^ ®s] 

n=3 

Therefore we rewrite equation (28) as 

(28a) a^^^T^(l) ^ag^^^^^l) + ^3 Tgd) + a^^ ^ ^ x^d) ^a^^^^x^d) + 

+ x^(2) + x^(2) + X 3 ( 2 ) + x^(2) + Xg(2) > B^ 
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Thus the five restraint equations generated from the second 
year requirements are equations (24a) , (25a) , (26a) , (27a) , and (28a) . These 
equations, in addition to the five equations (lla) through (lie) , provide the 
ten linear restraint equations for the two year models They serve to provide 
the inputs (A, B, cY ) to the linear program which will minimize the training 
input during both years of the two year projection, while meeting requirements 
both years . 

These inputs are: 



1 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


0 


1 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


^,1 


^8,2 


^9,3 


^10,4 


^11,5 


1 


1 


1 


1 


1 


^8,1 


^8,2 


^9,3 


^10,4 


^11,5 


0 


1 


1 


1 


1 


0 


^9,2 


^9,3 


^10,4 


^11,5 


0 


0 


1 


1 


1 


0 


0 


®10,3 


^10,4 


^11,5 


0 


0 


0 


1 


1 


0 


0 


0 


^11,4 


^11,5 


0 


0 


0 


0 


1 



Note; The rows of A (i = 1 ,2 , . „ . , 10 ) represent, in order, the ten 
restraint equations (lla-lle, 24a , 25a ,26a , 27a , 28a ) , The 
elements of each row are the coefficients of the ten unknown 
training Inputs Tj^(ni) (k = 1 , 2 , . . . , 5; m = 1 , 2 ) . 
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10 



where are defined by equations found in the derivations of the restraint 
equations; 



c=[cj = (Oj, O^, C 3 , C^, Cj) 



where the are defined in Chapter II; 

Y=[T^(m)] m=l,2; k=l,2,..„,5 

where the Tj^(m) represent the ten training inputs during the two years, 
and the objective function is 

2 5 

minimize cY = min ^ ^ c^ 

n^l k=l ^ 



With the ten training inputs having been determined by 
the linear program, equations ( 22 a) through ( 22 f) are used to determine the 
training inputs to each of the 72 cells during both years of projection, 

T (m), ( m = 1 ,2; n = 1 ,2 , . . „ ,72) . Knowing the values of all T (m) , the 
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1 





projection equation (23) is then used to detemine each P (2) , (n = 1,2 , . . .,72) , 

n 

and hence the total projected inventory after two years, P (2) . 
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b. Five Year Model 



This section of Appendix B lists the explicit expressions 
for the inputs to the linear program . 

1 . Matrix of Coefficients , A = (a ,^) 

The matrix of coefficients, A = (a_) , for the five year 
model is shown in Figure B-1 . 

The computational formulae for describing the matrix 
A = (a_) for the five year model are as follows; 

L^n]^ ^8.1 

n=3 

24 29 

■3.1- ■''•LL”. ' • *1C,3 

n=21 n— 27 

36 42 

42 48 

"lO,4=l/'[LUn ^ "ll,4 

n=39 n=45 

54 60 

^ 1 . 4 = ^Ijn] 

^2,1= 1/2[M[^7,i'"8,i] " "l3,l 
n=4 

^3,1= ?;[n]L"7.r^8.rj ^^4,1 

n=10 n=10 

"13,2= ‘4?Un^?Un]K.2-^9,2]^" 
n=ll n=l6 
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14,2 





i 

I 




\ 




I 



14 , r 

n=22 



® 14 , 2 " 

^ 14 , 3 " 
^ 15 , 3 " 
^ 15 , 4 " 
^ 16 , 4 " 
^16,5^ 
^17, 1" 

^18, 1" 
^18,2^ 
^19, 1" 
^19,2" 
^19,3^ 



:|jn]"9,2 " l/6[|w . I wj[a3_2- ] 

25 30 

n— 22 n— 2o 

^?jn]"l0,3^ " 10 , 3 ] 

n=34 n=40 n=34 n=40 

43 48 

‘^■LLvt»n]-n,i 

n=52 n^8 

I}^n][^l2,r ^13,l] ^18,1 

n=5 

f ^n]"l3 , 1^ ‘/'L 1^] 1 12 . r ^3 , 1 ] ^ ^9 , 1 

n=ll n^l 

l/=[^V?,%]["l3,2-^4,2]^ ^9,2 
n=l2 n=l/ 

‘/•iX'J •..,.* ■/<L|.-'Jt-.,..-»,.l 

26 30 25 30 



20,2 



a 



38 43 

20,2=l4Ijn"L^] 

n=36 n=4l 
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20,3 

^20,4^ 

^21,4" 

^21, s" 

^22, l" 
^23, l" 

^23, 2^' 
^24, l" 
^24,2" 
^24,3" 
^25, r 

^25,2^ 

^25,3" 

^25,4" 



38 44 37 43 

l/aUJntL^>i5,3^ l37n;Z^Wn]K4,3- ^15,3. 



n=35 n=4l 
_ 44 48 



15,4 16,4 



56 62 56 61 

^6 [^17, r ^18, l] ^23,1 

l/^[X''"Jtl 7 .r"l 8 ,l]" ^24.1 
Z '^ 18 ][^ 18 , 2 ‘ *19,2] ® 24,2 

n=13 

!/■*[ L^n] ^9 . 1^ Z,%] 1 18 , r ^9 , 1 ] ^ ^25 , 1 

n=24 n=24 

Ly '' 30 ] ^19,2" 1/4 Ly ’^3o]h8.2-^9.2]^ 

n=25 n=25 

27 

^3o]lf 19 , 3 " ^ 20 , 3 ] ^25,3 

39 

^9,1 

n=36 

39 48 39 49 

l/5f y Y„+y Y^la__ _+ 1/6F y W + V W^Ta^- o" ^on 9 

^n 437 %^ 47 ^^ ^n 4 s 7%^47 

• 45 38 43 

Z,y . l/7[ y y y wj[a -a^ ] 



‘n=35 n=4l 

Lyn"^ ^48t20,4“ ^ 21 , 4 ] ‘^^26,4 



■n=42 
. 45 



^25,2 



] 

^26,3 
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I 



20,3 



26 



.3= 



n=49 



57 63 57 y 

,4= l/lLyI^]^21,4- ^^il^yljn]b20.r^nj 



n=54 n=60 

57 63 

^26,5= 1/®[I U„-Zjn]^21,J 
n=54 n=60 



2„ Restraint Vector, B = (B^) 

The restraint vector for the five year model has 25 
elements, five for each year. The elements are computed for each year 
by a method similar to that used in the one and two year models . 

Bo 



Second Year 



12 



16 

n. 

-\7 



21 

^22 



2 6 



Third Year 



Fourth Year 



Fifth Year 
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3 Cost Vector, ~ 



The cost vector is defined in Chapter II. 

4. Unknown Training Input, Y 

Y = Tj^(m) J , m = 1 , 2 , . , . , 5; k = 1 , 2 , _ . , 5; 
where represents the 25 training inputs during the 

five years , 

5. Objective Function, min cY 



min cY = min 



1 , 1 , 



m=I k=I 
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APPENDIX C 



SENSITIVITY ANALYSIS 

A cor^iputer model like PIPE is as effective and accurate 
as the input parameters that it uses. The accuracy required of 
the input parameters, however, may be dictated by the effect that 
changes in these input parameters have on the resultant solution. 
Sensitivity analysis is the method by which this effect is evaluated. 

PIPE FIVE was evaluated for sensitivity by successively 
operating the model with all parameters at "standard" values, 

(as explained in Section II C), except one. This one was, in turn, 
attrition rate, promotion rate, reenlistment rate, cost coefficients, 
and control factor. The effects of varying these parameters, one 
at a time, is detailed in the following sections. 

1. Variation of Attrition Rate 

The five-year model was operated with all parameters at 
standard values except attrition rate, which was varied from 0. 003 
to 0. 027 in increments of 0. 002. Figure C-1 illustrates the results 
of these projections, the fpllowing of which are significant: 

a. Total five-year cost increases linearly as attrition 
rate increases. This is intuitively the expected result, since the 
program should have higher cost if more men are leaving the program 
through attrition. 
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b. The curve for the total five-year input has a sharp 
break in the interval where the attrition rate is 0. 013 to 0. 015. This 
demonstrates the fact that for the minimum cost solution to the 
problem, the inputs of personnel consist of men in the lower (less 
expensive) pay grades. Near an attrition of 0. 015, however, it 
becomes less expensive to put in less men in the higher pay grades, 
(see Figure C-2A). This "breakpoint" is not intuitively obvious to 
the personnel planner, but is clearly delineated by PIPE. 

C. Attrition rate has an effect on the number of men 
required by the program during any given year of projection. 

Figure C-2B demonstrates that for attrition rates of 0. 003 to 
0. 013, about twice as many men are put into the program during 
the third year as in the second and fourth years. For attrition 
rates of 0. 015 to 0. 027, however, the fourth year of projection 
requires double the input of men compared to the second and third 
years. 
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FIGURE C-l 



127 



FIVE YEAR TOTAL COST 
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ATTRITION RATE 
FIGURE C-2 



2 . 



Variation of Promotion Rate 



a. Figure C-3 illustrates that as the promotion rate is 
increased for all pay grades, the total cost of meeting requirements 
decreases while the total training input fluctuates. The reasons for 
this fluctuation are illustrated in Figure C-4. 

b. Figure C-4 illustrates that as the promotion rate increases, 
the training inputs to E-4, E-5, and E-6 decrease while that of 

E-3 increases. This happens because the higher promotion rate 
enables the requirements for the higher pay grade to b e met by 
the increasing portion of the E-3 input that is promoted during the 
year of projection. 

c. The training input solution for pay grade E-7 is zero for 
all promotion rates considered. 

d. It appears from Figure C-3 that the promotion scheme 
representing the minimum total training input is about at the values 
of promotion rate used as "standard", (i. e. , at delta gamma equal 
to zero). Thus it may be the case that, for minimum input of 
personnel over the five-year period, promotion rate should be at 
the "standard" value, but for minimum cost, promotion rate should 
be as high as possible. 
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FIVE YEAR TOTAL TRAINING INPUT 



CHANGE IN PROMOTION RATE 
VS 

FIVE YEAR TOTAL TRAINING INPUT 
AND 

FIVE YEAR TOTAL COST 




FIGURE C-3 
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CHANGE IN PROMOTION RATE 
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CHANGE IN PROMOTION RATE FROM STANDARD VALUE, (AT) 
FIGURE C-4 



3. 



Variation of Reenlistment Rate 



Sensitivity of the model to changes in reenlistment rate at the 
end of the first enlistment was performed on all pay grades since it 
was felt that this was the most critical time in the reenlistment 
problem. Reenlistment rate was varied in increments of 0. 05 to 
values of -0. 30 and +0. 30 from the standard value. In the interval 
from the standard value to +0. 12, a fluctuation in the total training 
input was noted. Further analysis was therefore made of this interval, 
with reenlistment rate changed by increments of 0. 005 in this region. 

a. Figure C-5 illustrates that as reenlistment rate is in- 
creased for all pay grades, the total cost of meeting requirements 
decreases while the total training input fluctuates. The reasons for 
these fluctuations are illustrated in Figure C-6. 

b. Figure C-6 illustrates that as the reenlistment rate in- 
creases, the training inputs to pay grades E-4 and E-5 decrease while 
that of E-3 increases. This happens because the higher reenlistment 
rate enables the requirements for E-4 and E-5 to be met by the in- 
creasing portion of the E-3 input that is retained during the years of 
projection. 

c. The training input solution for pay grades E-6 and E-7 
are zero for all reenlistment rates considered. 

d. Figure C-5 illustrates that the minimum cost solution 
and the minimum training input solution are both found at the point 
of maximum reenlistment rate. 
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CHANGE IN REENLISTMENT RATE FROM STANDARD (A8) 



FIGURE C-5 
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FIVE-YEAR TOTAL COST 




CHANGE IN REENLISTMENT RATE 
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hA8 = -0.30 to +0.30 }-A8=-0.30to+0.30-H kA85-0.30to+0.30-H 

CHANGE OF REENLISTMENT RATE FROM STANDARD (A8) 



4. Cost Coefficients Changed 

The cost coefficients in PIPE FIVE represent the relative 
costs between the different pay grades considered. These cost 
coefficients were varied by changing the relative cost or value of 
one rate as compared to another. Thus an increase in relative 
costs indicates that one rate has become comparatively more ex- 
pensive to use when compared to a lower rate than had previously 
been the case. 

a. The total five-year cost of the program increases almost 
linearly as relative cost increases, (Figure C-7). 

b. Figure C-7 shows that the total five-year training input 
increases linearly with increase in relative costs to a certain point, 
(about 0. 15), when a plateau is reached. Beyond a relative cost 
within pay grades of 15%, the cost of the higher rates becomes so 
prohibitive that the minimum cost solution dictates that all input be 

in the pay grades E-3, E-4 and a smaller amount of E-5, (Figure C-8). 
At this point, since a certain number of personnel are needed to meet 
requirements, the same numbers of "cheap" personnel are brought 
into the program. Thus a further increase of relative cost would 
have no effect on the slope of the training input line. 
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PROPORTIONAL INCREASE IN COST PER PAY GRADE 



5. 



Variation of Control Factor 



The Control Factor of PIPE FIVE is a restraint on the 
number of personnel that the model may carry in excess of 
requirements during any one year in order to derive the minimum 
cost solution over a five year period of projection, (see 
Appendix D). Thus a control factor of 0. 0 requires the model 
to exactly meet requirements each year, whereas a control 
factor of 10% allows as much as 10% excess over requirements 
in any pay grade in any year as long as requirements are met. 
Control factors from 0. 0 to 50. 0% were used to test the 
sensitivity of the model to this parameter. 

a. Figure C-9 demonstrates that the total training 
input does not vary as the control factor is increased and more 
excesses are allowed. This is because minimum requirements 
do not change, and must be met without regard to the control 
factor. 

b. Figure C-9 illustrates that the total five-year cost 
decreases as the ceiling on excesses is relaxed. The reasons 
for this may be seen in Figure C-10. As more excesses are 
allowed in the training inputs, more E-3's are procured 
because they are cheaper than the E-4's, whose number decreases. 
The input of E-5's remains constant. 
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c. Figure C-11 shows the effect that an increase in 
control factor has on the year of input of the personnel. As more 
excesses are allowed, the number of people who enter in the 
fourth year of projection increases, while the number of those 
who enter in during the fifth year decreases. 
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FIVE-YEAR TOTAL TRAINING INPUT 
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APPENDIX D 



USER MANUAL 

1. INTRODUCTION AND GENERAL DESCRIPTION 

This Appendix is designed to provide the personnel planner 
with the information necessary to use the computer programs without 
a detailed knowledge of either the mathematical model from which 
they were derived or the techniques used in their construction. 

The computer models described in this Appendix are the One 
Year Model (PIPE ONE), and the Five Year Model (PIPE FIVE). 

The One Year Model is essentially a rnechanization of the techniques 
now used by personnel planners in SUPERS to project inventories 
using adding machines and desk calculators. The Five Year Model 
performs this same function but with two very important sophisti- 
cations; a computational technique is used which: 

a. allows the model to "look ahead" at future requirements 
and the necessary training inputs that these requirements 
demand; 

b. adds these training inputs to the inventory over the five 
year period in such a manner that the total cost of these 
inputs is minimized while meeting requirements each year. 

The ability of the computer to account for the future and to attach the 
notion of "cost" to the computed results are the significant additions 
to the planners' "bag of tools" offered by the Five Year Model. 
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The solution generated by a computer to a problem as complex 
as the planning of future enlisted inventories should not be regarded 
as a set of final precise figures, but rather as a guide to decision. 

The computer models are an attempt to simulate actual conditions 
(i. e. , the natural environment) to a degree that will yield predictions 
accurate enough to make the application of the model worthwhile. 

An objective in formulating these models was that they would 
produce a solution of accuracy comparable to the accuracy of the 
input data, and provide a rapid analytical means of solving enlisted 
projection problems using data familiar to the personnel planner. 

2. INPUTS TO THE MODEL 

This section discusses the inputs to the computer models which 
must be determined by the user. After the data has been assembled 
in the manner described in the following paragraphs it must be trans- 
ferred to data sheets. Sample data sheets filled in with representa- 
tive data for each computer input are provided. These data sheets 
are designed to assist the planner in preparing the data in a manner 
in which it can be readily interpreted by the computer key punch 
operators. 
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The inputs to the computer program are: 

a. Initial Inventory 

b. Requirements 

C. Training Input Distribution 

d. Cost Coefficients (only PIPE FIVE) 

e. Attrition Rate 

f. Retirement Rate 

g. Reenlistment Rate 

h. Promotion Rate 

i. Control Factor (only PIPE FIVE) 

j. Serial Number 

These inputs are examined in detail in the following sections. 

a. Initial Inventory: These models will project enlisted 

inventories of personnel with the following characteristics: 

(1) similar rates of attrition, reenlistment, retire- 
ment and promotion; 

(2) similar requirements to be met 

For example, machinist mates and electronic technicians could 
conceivably satisfy the first criterion because of similar character- 
istics of promotion, attrition, reenlistment and retirement; however, 
they can not be projected as one group because they obviously do not 
fill similar billets. 
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Figure II- 3 
Inventory Matrix 



In order to be used by the computer the initial inventory must 
be subdivided by pay grade, obligated service, and current en- 
listment into a matrix of 72 cells as shown in Figure II-3, repeated 
here for the user's convenience. The six rows of the matrix represent 
years of obligated service and the twelve columns are divided among 
five pay grades. For the purpose of this matrix, E-1, E-2 and E-3 
are compressed into one group, E-3. Similarly, E-7, E-8, and E-9 
are considered together as E-7. Each column represents an enlist- 
ment, first, second, third, fourth or fifth. Note that there is pro- 
vision for only one enlistment for E-3, two enlistments for E-4 and 
E-5, and three enlistments for E-6. The assumption was made that 
for critical ratings the number of personnel who did not fall within 
the limits of these pay grades would be inconsequential and could be 
neglected. Note also that there is no provision for E-7 in their first 
enlistment. The assumption is made that the numbers of people in 
this category (first enlistment chiefs) could be safely ignored. 

The criterion for allocating individuals of the inventory into 
the appropriate rows and columns (i. e. , cells) of the inventory 
matrix are given below: 

a. Row Criterion 

The row position of the members of the inventory is 
determined by the number of years of obligated service of an 
individual. Obligated service for this purpose is defined as the 
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difference between the calendar year of an individual's expiration 
of active obligated service and the present calendar year, plus one. 

Let: n = (calendar year of expiration active obligated service) 

- (present calendar year) 

Then: row position = n + 1 

Example: a person in a calendar year 1964 inventory whose expiration 
of active service occurs in calendar year 1967 would be placed in row 
number four. That is, 

n = 1967 - 1964 = 3, so that 
row position = n+ l= 3 + 1 = 4 
(2) Column Criterion 

The column position of the members of the inventory 
is determined by their pay grade and current enlistment. 

The following table summarizes the column breakdown: 



COLUMN 


PAY GRADE 


CURRENT 

ENLISTMENT 


LENGTH OF SERVICE 
IN YEARS 


1 


E-3''^ 


First 


0 to 6 


2 


E-4 


First 


0 to 6 


3 


E-4 


Second 


7 to 12 


4 


E-5 


First 


0 to 6 


5 


E-5 


Second 


7 to 12 


6 


E-6 


First 


0 to 6 


7 


E-6 


Second 


7 to 12 


8 


E-6 


Third 


13 to 18 


9 


E-7>!^>!' 


Second 


7 to 12 


10 


E-7'!'^« 


Third 


13 to 18 


11 


E-7^'^' 


Fourth 


19 to 24 


12 


E-7^1^^'* 


Fifth 


25 to 30 



Notes: * includes pay grades E-1, 2, and 3 

includes pay grades E-7, 8, and 9 
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After the inventory has been assembled in the manner described 
above, it must be transferred to data sheets. Figure D-1, Initial 
Inventory, is the data sheet used. The numbers above each box cor- 
respond to the column numbers on a standard 80 -column IBM punch 
card onto which the item of data is punched. These numbers do not 
concern the planner. The appropriate data is placed in the boxes, 
one digit per box. There is one precaution in filling out the data 
sheets; whole numbers without decimal points must be right adjusted; 
percentages and decimal fractions are left adjusted. That is to say, 
a particular group of data boxes must be filled so that there are no 
blanks on the right or left as the case may be. A blank box will be 
interpreted as a zero by the computer. In both right and left 
adjusted data the user may indicate the proper location of the decimal 
by placing a decimal point in the appropriate box, keeping in mind 
that blanks are read as zeros. 

A completed Initial Inventory Data Sheet, Figure D-1, represents 
one column of the inventory matrix. Twelve of these data sheets are 
required for the initial inventory, each sheet representing a column 
of the inventory matrix. The data are placed in the boxes as shown 
in the figure. 
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Column No. ^ INITIAL INVENTORY - INPUT DATA SHEET Rate 
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b. Requirements: As used by the computer programs re- 
quirements must: 

(1) be expressed by pay grade and, 

(2) correspond to the initial inventory. 

These requirements, by pay grade, must be defined for each year 
of projection, and may differ from year to year. 

A completed Requirements Data Sheet for a five year projection 
is found in Figure D-2. One completed data sheet is required for 
a five year projection. 

c. Training Input Distribution: As stated in the introduction to 
this Appendix, one of the outputs of the computer program is the 
"training inputs" which are required to maintain the projected 
inventory at a strength necessary to meet requirements. The 
distribution of the training input in the inventory matrix, i. e. , to 
which cells this training input will actually be added is determined 
by the personnel planner. He must identify the amount of obligated 
service and the enlistment number for the computed training input to 
each pay grade. The model is constructed so that input may be made 
to any cell which represents personnel with less than 18 years of 
service, (all columns except 11 and 12). The user must determine 
for each pay grade and for each column within each pay grade: 

(1) If there will be a training input in a particular 
column and. 
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REQUIREMENTS - INPUT DATA SHEET Rate 
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FIGURE D-2 



(2) If there is an input into which particular cells of the 
column the input will be allowed. 

Once the training input distribution has been established the computer 
solves the problem for the number of people who must be added to a 
particular pay grade to maintain requirements, and then divides 
this input equally among the cells of the pay grade that have been 
designated by the user. 

Looking at the data input sheet, Figure D-3, for the Training 
Input Distribution, the terms "upper" and "lower" limits of a cell 
must be defined. The "upper" limit is the largest cell number in 
the column into which training input will be added, and the "lower" 
limit is the smallest cell number. For example, if in column #3 
the training input is to be placed in cells 13, 14, 15 and 16, then the 
upper limit would be 16 and the lower limit 13. In determining the 
limits for any particular column there is the restriction that in any 
one coliomn the cells into which training input may be assigned must 
be consecutive. It is important that each limit be filled either with 
a number corresponding to a cell number or a zero, which indicates 
no input. If input is to be allowed in only one cell of a particular 
column then this cell number must appear as the "lower" limit and 
a zero must be entered as the "upper" limit. 
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TRAINING INPUT DISTRIBUTION - INPUT DATA SHEET Rate 



cc! - 
W ^ 
P. 5 
a, s 






W H 

o a 
j a 






pcj 

3 H 
a 5 

a a 
D a 

Pi 

w H 

la 

a a 



o a 
u a 



nO 



v> 



|n 1 00 I cn| 00 







^ w 


















a a 




^ u 




^ a 




D 






a pej ► 






a > 


a ^ 






0 K < 


D a 


( 


D M 


0 a 






U H ( 


j a 


U w 


U H 


pel - 




















W ^ 




















a 5 so 

a a 




vO 


- 


vO 


(O 


vD 


0 


sO 


pn 


^ un 




in 




le 


rvl 


m 




Lfl 


v> 


ttj 




















a H fT^ 




CO 


(T 


CO 


O 


CO 


O 


CO 


O 


o ^ 

“S'" 




(M 




oo 








(M 


Vrt 




o 



O a; 

5 | 

C (U 
■■^ X 

QJ 

c a 



>H < 

< a 
a o 



154 



FIGURE D-3 




The following examples indicate how these cell limits are 
determined and also suggest various training input distributions 
that might be used: 

Example (1): Training input is to be allowed in all grades, but the 

input must have four or more years of obligated service. The limits 
for the respective columns are: 



Column Number Lower Limit 


Upper Limit 


1 


1 


3 


2 


7 


9 


3 


13 


15 


4 


19 


21 


5 


25 


27 


6 


31 


33 


7 


37 


39 


8 


43 


45 


9 


49 


51 


10 


55 


57 


Example (2): Training input 


is to be allowed in 


the first enlistment 


only and the input must have 


more than four years of obligated 


service. The limits for the 


respective columns 


are: 


Column Number Lower Limit 


Upper Limit 


1 


1 


2 


2 


7 


8 


3 


0 


0 


4 


19 


20 


5 


0 


0 


6 


31 


32 


7 


0 


0 


8 


0 


0 


9 


0 


0 


10 


0 


0 
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Example (3); Training input is to be allowed in the first enlistment 
and the input must have more than four years of obligated service 
except that an E-6 with 6 years of obligated service may be accepted. 
The limits of the respective columns are the same as in Example (2) 
above except for column #7 where the lower limit would be 37 and 
the upper limit would be 0. 

Example (4): Training input is to be allowed for personnel with 

between 4 and 8 years of service. The limits for the respective 
columns are: 



Column Number 


Lower Limit 


Upper Limit 


1 


4 


6 


2 


10 


12 


3 


13 


14 


4 


22 


24 


5 


25 


26 


6 


34 


36 


7 


37 


38 


8 


0 


0 


9 


49 


50 


10 


0 


0 



Example (5): No training input will be allowed for E-6 and E-7. 

(This might be a situation where no matter what solution the computer 
might generate which says "add E-6 and E-7" these personnel may be 
unavailable and hence such a solution is unacceptable. ) If the train- 
ing input to a particular grade is suppressed the user must accept 
the possibility that the numerical requirements for this particular 
grade may not be met; however, overall numerical requirement for 
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the inventory will always be met, provided that the input to the 
lowest grade E-3 is not that pay grade that is suppressed. In this 
example, if the requirements for E-6 and E-7 are not met they will 
be numerically filled by E-5 and lower pay grades. For this 
example, the "Lower" and "Upper" limits for columns #6 through 
#10 would be zero. 

d. Cost Coefficients: The cost coefficients provide a means 

whereby the computer model can consider the relative costs of ad- 
ding personnel of different pay grades to the inventory. The cost 
coefficients measure the relative difficulty of obtaining personnel 
of one pay grade as opposed to another. They may be expressed 
in terms of dollars, or as ratios which represent procurement 
difficulties. J£ expressed in terms of dollars, the planner may 
consider training costs, pay and allowances, etc. In the process 
of solving for the training input to meet requirements, the computer 
considers the cost of adding this input to the various pay grades and 
then solves for the least expensive solution, based on the cost co- 
efficients determined by the user. 

Five cost coefficients are required, one for each pay grade. 
Figure D-4 is a sample cost coefficient data sheet completed with 
the cost coefficients used by the authors in their evaluation of the 
Five Year Model. Note that these coefficients may use the decimal 
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COST COEFFICIENTS, CONTROL FACTOR, AND SERIAL NUMBER 
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point. If the data is expressed in decimal form, as in the example, 
Figure D-4, the user must place the decimal point in the proper box. 

If no decimal is used, the numbers must be right adjusted as 
explained previously. 

e. Attrition Rate: Associated with the passage of each in- 
crement of time is the likelihood that a certain fraction of the in- 
ventory will be lost or "attrited" due to unforseen circumstances. 

In order to account for this "unplanned" loss of personnel, an 
attrition rate is applied to the inventory in the model. This rate is 
composed of the summation of specific types of losses expressed as 
a percentage reduction to be applied to the inventory each year. 

Attrition rate consists of the following types of personnel losses: 
(The numbers in parentheses following the various types of losses 
are the specific item numbers for attrition of active duty enlisted 
personnel as defined in NAVPERS 15658, Navy and Marine Corps 
Military Personnel Statistics . ) 

(1) Death (#1290) 

(2) Physical Disability (#1120) 

(3) Discharge by Undesirable Discharge (#1000) 

(4) Discharge by Bad Conduct Discharge (#1040) 

(5) Discharge by Dishonorable Discharge (#1050) 

(6) Discharge by Reason of Dependency or Hardship (#890) 

(7) Intra-Navy Transfer (#1250) 

(8) Miscellaneous 

The attrition rate data sheet. Figure D-5, shows that an 
attrition rate is entered for each cell of the matrix. The attrition 
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Column No. 7 ATTRITION RATE - INPUT DATA SHEET Rate 
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FIGURE D- 



rate is a decimal fraction which is always less than one. In entering 
the data the user may either indicate the position of the decimal point 
or the data must be left adjusted so that it is properly interpreted 
by the computer. 

f. Retirement Rate: The retirement rate is defined as the 

percentage of personnel in the inventory who are lost through 
voluntary retirement by reason of having completed 20 or more 
years of active service. Since only 11 cells of the inventory, (in 
pay grade E-7, cells 62 through 72), represent a total length of 
service of 20 years or more, only those cells have an associated 
retirement rate. A retirement rate must be assigned to each of 
these 11 cells. For example, if fifty percent of the E-7's who start 
the twentieth year are expected to retire within one year then the 
retirement rate for 62 would be 0. 5 (point five). A completed 
sample retirement data sheet is shown in Figure D-6. 

g. Reenlistment Rate: Reenlistment rate is the ratio of the 
number of men reenlisting to those eligible to reenlist, expressed 

as a percentage. The reenlistment rate is applied at the termination 
of a person's obligated service. Looking at the inventory matrix 
in Figure II-3, the top row represents personnel with one year of 
obligated service. In each year of projection the personnel in this 
row are subject to a reenlistment rate. The percentage who reenlist 
are added to the inventory of the next higher enlistment in the cell 
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RETIREMENT RATE - INPUT DATA SHEET Rate 
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FIGURE D-6 



which represents 6 years of obligated service; the percentage that 
fail to reenlist are lost from the inventory. For personnel in any 
cell except those in the top row of the inventory matrix the value 
of the reenlistment rate is 100%. 

The planner must determine twelve reenlistment rates; E-3 
at the end of the first enlistment, E-4 and E-5 at the end of the 
first and second enlistments, E-6 at the end of the first, second, 
and third enlistments, and E-7 at the end of the second, third, 
fourth, and fifth enlistments. One completed data sheet is required 
for each projection as per the sample. Figure D-7. 

1. Promotion Rate: Promotion rate is defined as the per- 
centage of men promoted from one enlisted pay grade to the next 
higher enlisted pay grade in one year. Promotion is limited to 
intra-enlisted promotions, and since the inventory matrix considers 
E-7, E-8 and E-9 as one group, (E-7), only cells in pay grades E-6 
and junior have an associated promotion rate. Furthermore, pay 
grades E-1, E-2 and E-3 are compressed into a single pay grade, 
E-3, so that the lowest promotions are to E-4. Promotion rates 
for the individual cells may be determined from statistical data of 
past promotions, from current promotion plans or from any other 
similar source. 

Except for cells 31 through 3 5, all cells below cell number 62 
are assigned a promotion rate between and including zero to 100%. 
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REENLISTMENT RATE - INPUT DATA SHEET Rate 
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FIGURE D-' 



In general, promotion is between cells in the same enlistment, i. e. , 
an E-4 in the first enlistment is promoted to a cell in pay grade E-5 
which is in the first enlistment column of that pay grade. (The only 
exception is at the reenlistment point where a may may be considered 
to be promoted and reenlisted at the same time in which case he 
would move to the next higher pay grade in the next higher enlist- 
ment column. ) Cells 31 to 35 which represent E-6 in their first 
enlistment have zero rates of promotion because pay grade E-7 in 
the model starts in the second enlistment; hence movement from the 
first enlistment of E-6 is impossible. 

Figure D-8 is a completed sample Promotion Rate data sheet. 
Eight data sheets are required, one for each column. 

i. Control Factor: As explained in paragraph d. , "Cost 
Coefficients, " the computer solution for the training input to meet 
requirements is a minimum cost solution. This means that the 
computer has solved the problem without regard to how many people 
are added in any one year except for the restriction that the numerical 
requirements must be satisfied. Therefore, unless the computer is 
somehow restrained it will not consider the actual yearly distri- 
bution of the input, nor will there be any limit on the excess of 
requirements that may occur in any one year of the projection. 
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Column No. ^ PROMOTION RATE - INPUT DATA SHEET Rate 
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The minimum cost solution for the overall period is simply 
the least expensive means of meeting requirements, through train- 
ing input, without regard to the resultant distribution of the inventory 
in any particular year. As an example, the minimum cost solution 
for the five year period may call for all the total input to be 
affected in the first two years. Conceptually this answer may be 
correct, as this may be the least expensive solution. However, as 
a practical procedure this solution may be unworkable because an 
even flow over the five year span may be required by the training 
facility capacity. Furthermore, the requirements in this example 
would probably be greatly exceeded in the first two years, tapering 
off to just meeting requirements in the fifth year. 

To limit the number of personnel in any one year of the projection 
that may be in excess of requirements, an input called the "control 
factor" is provided. This is a number which controls the amount by 
which requirements may be exceeded in solving for a minimum cost 
solution. The control factor is defined as the percentage of the total 
requirements by which the inventory may exceed requirements in 
any one year. For example, if the control factor is 10% and total 
requirements are for 800 personnel (i. e. , in all pay grades), the 
solution is limited to a total of 880 people in the inventory even if 
a more economical solution would dictate a greater number in a 



particular year. 



The smaller the control factor, the more restrictive and 



expensive will be the solution. The most restrictive value of the 
control factor is zero, which causes the computer to fill requirements 
as vacancies occur. This is the maximum cost solution in that the 
ability of the model to anticipate future losses is not utilized and 
the advantage of flexibility of solution has been removed. It was 
determined through analysis that using a control factor of 100% 
effectively removes any restriction on the amount that inputs may 
exceed requirements. In using the five year model it is recommended 
that a control factor of 100% be used to first find the minimum cost 
solution. If this solution is not satisfactory with respect to yearly 
input the control factor can be adjusted so that the requirements in 
any one year will not be excessive. The control factor input is 
illustrated in Figure D-4. 

j. Serial Number: This input provides the user with a means 

of serializing a computer run. This number may take any value 
from 1 to 9999 but may not be a decimal. This input is entered on 
the data sheet for "Cost Coefficients", Figure D-4 
3. OUTPUTS OF THE MODEL 

The basic outputs of the model are: 

a. Projected inventory; 

b. Training input; 

c. Training input cost, (PIPE FIVE only). 



The computer is programmed to display these outputs in two forms. 
Yearly Inventory Summary, Figure D-9 and Five Year Summary, 
Figure D-10. These two output sheets are summarized in the fol- 
lowing paragraphs. 

a. Yearly Summary, Figure D-9. 

(1) In the left upper corner "RUN NO. 1", the numerical 
value is the run "index" or serial number which is 

a result of input number 10. 

(2) "INVENTORY PROJECTION YEAR NO. " is the 

year of projection that is represented by the in- 
ventory matrix. "Year No. 1" would represent 

the inventory at the end of the first projection year. 

(3) "TOTAL PROJECTED INVENTORY" represents 
the pay grade totals and grand total of the inventory 
matrix. 

(4) "TOTAL TRAINING INPUTS" summarizes the 
training inputs by pay grade and indicates an over- 
all total. 

(5) "EXCESS (SHORTAGE) OF REQ" indicates how the 
total projected inventory for a pay grade compares 
with requirements for that grade. A positive 
number indicates an excess of requirements, where 
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RUN NO. 1 INVENTORY PROJECTION YEAR 
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YEARLY SUMMARY 
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the projected total does not meet requirements, 
a minus number indicates the shortage. Only under 
most unusual circumstances will the actual inventory 
precisely fit the requirements structure. It is re- 
called that the computer is programmed to meet 
total requirements from the "top" down so that a 
shortage in E-3, for example, is made up for by 
an excess in a higher pay grade. 

(6) "TRAINING INPUT PER CELL" shows the training 
input into the inventory matrix by cell. 

b. Five Year Summary, Figure D-10, This Output Form 
repeats the EXCESS (SHORTAGE) OF REQ, TOTAL PROJECTED 
INVENTORY, and TOTAL TRAINING INPUTS totals from the one 
year summaries. The TOTAL FIVE YEAR TRAINING INPUT is the 
summation of the yearly training totals. The FIVE YEAR COST is 
the total cost of the training input. It is computed by multiplying 
the total five year training input for a pay grade by the cost co- 
efficient for that pay grade, and summing over all pay grades. 

4. SPEED OF COMPUTATION (CDC 1604 

a, PIPE ONE; 50 seconds to compile, then one second for 
each year of projection. 

b. PIPE FIVE; 2 minutes 50 seconds to compile, then 15 
seconds for each year of projection. 
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It is obvious that one of the major advantages of these computer 
programs is their speed, A planner faced with selecting a promotion 
rate for a particular grade can evaluate many plans and their affects 
in a very short period of time. The most efficient use of the com- 
puter, because of the relatively long time required by the compiling 
processes, is to make as many projections as possible at one time. 

5. USES FOR PIPE ONE 

a. PIPE ONE projects inventories for one year, computes 
the training input to meet total requirements (filling senior require- 
ments first), and may be cycled for any number of years to give a 
long range projection. The advantages of this program are that it 
takes little computer time and can be used to make projections 
beyond five years. The disadvantages when compared to PIPE FIVE 
are that it does not give the "cost" information and it merely fills 
vacancies as they occur. 

b. PIPE ONE, with a minor program change, can be made 
to project inventories with a constant training input into any cell or 
group of cells of the inventory matrix. In this manner proposed 
training input policies may be evaluated. 

c. The effects of various promotion and retirement policies 
on the projected inventories and training inputs may be evaluated. 
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d. The effects of an increase (or decrease) in reenlistment 
rates may be evaluated. 

e. PIPE ONE, with a minor program change, can be used 
to build an enlisted "model structure" for any number of years for 
any particular set of promotion, attrition, reenlistment and retire- 
ment rates. This may be accomplished by a program change which 
would provide a training input of 100 people in cell #1 in each year 
of projection. Initial inventory and requirements are set at zero. 

The projected inventory for the last year of projection will show for 
each cell the percentage of the initial inventory, i. e. , 100 people 
still on board after the passage of time. 

6. USES FOR PIPE FIVE 

a. PIPE FIVE projects inventories for five years, computes 
the training input so as to meet total requirements in the least ex- 
pensive manner and determines the cost of this training input. This 
program takes into account future losses and changes in requirements 
in computing the training input. The program anticiaptes future needs 
while meeting current requirements and hence may compute a train- 
ing input that causes the total requirements to be exceeded in one or 
more of the years of projection. The amount that requirements may 
be exceeded in any year may be controlled by the user. 
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b. As in PIPE ONE, PIPE FIVE can evaluate the effects of 
promotion rate, reenlistment, and retirement rate on the inventory 
and on the training inputs required to insure that the inventory meets 
requirements. At the same time it can evaluate the changes in the 
cost of the training input for various changes in these rates. 

c. By adjusting the "control factor" the planner can examine 
a range of "cost" solutions and use the solution which is the best 
compromise between training cost, uniform yearly training input, 
and the number of people in excess of requirements. 
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COMPUTER PROGRAM LISTINGS 
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1 . 



PIPE ONE 



The following listing for PIPE ONE is written in FORTRAN 60. 
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SUBROUTINE FACTOR 01170 
ODIMENSICN X(72) , ALFA(72) ,BETA{72),GAMJ^A(72) ,DELTA(72) , 01180 
1P(5,72).T(5,72),TAU(5,5) ,REC(5,5),R{ 5,5) ,TEMP(5,6), 01190 
2C(5,72),D{5,7 2),CK(5,72),SP(5,6),ST( 5,6),EXR(5,6), 01200 
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68 TX(5) = 1./(TA(9) + TA( 10) ) OUUO 

69 CONTINUE 01450 

RETURN 01460 

END . - 01470 



SUBROUTINE EXCESS 01500 
OOIMENSION X{72),ALFA(72),BETA(72),GAMMM72),DELTA(72), 01510 
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ODIMENSION X( 72 ) ,ALFA( 72 ),BETA( 72 ),GAMf^A( 72 ) ,DELTA( 72 ), 01710 
1 P( 5 , 72 ),T( 5 , 72 ) ,TAU( 5 , 5 ) ,REC( 5 , 51 , R( 5 , 5 ) , TEMP ( 5 , 6 ) , 01720 
2 C( 5 , 72 ) ,D( 5 , 72 ),CK( 5 , 72 ) ,SP( 5 , 6 ),ST( 5 , 6 ) ,EXR( 5 , 6 ) , 01730 
3 TX( 5 ),TA( 10 ),L( 20 ) 01740 
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ODIMENSIGN X ( 7 2 ) , ALF A ( 72 ) , B ET A ( 72 ) , GAMF A ( 72 ) , DEL T A ( 72 ) , 02060 
1P(5,72»,T(5,72) , TAU ( 5 , 5 ) , R EG ( 5 , 5 ) , Rt 5 , 5 ) , TEMP ( 5 . 6 ) , 02070 
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ODIMENSION X(72) ,ALFA(T2),BETA(72),GAMMA(72),DELTA(72), 02290 
1P(5,72),T(5,72),TAU(5,5) ,REC(5,5),R( 5,5) ,TEMP(5,6), 02300 
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2 . 



PIPE FIVE 



The following listing for PIPE FIVE is written in FORTRAN 60. 
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GDI MENS ION X(72) , ALFA( 72 ), BETA ( 72 ), GAMMA! 72 ) , DEL T A ( 72 ) , UM B ( 72 ) , C0030 
1P(5,72),T(5,72) ,TAU(5,5) ,REC(5,5),R{ 5,5) ,TX< 5),TA{ 10) ,CC (5) , 00040 
2C(5.72),D(5,72),CK(5,72),SP{5,6),ST( 5,6) ,EXR(5,6). 00050 
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1300 FORMAT (14) 00300 

DO 1400 N=1 ,72 0031C 

1400 UMB(N)=X(N) 00320 
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COMPUTE CONSTANTS 



o o o o 

r- CVJ lO J- 
O' O* O' o* 

o o o o 
o o o o 



o o o o o ooooooooooooooooooooooooooo 
Lf>’or-coo'0*— c\j(v)ifmsot^coo'o*— csj'n-ifin'Oh-ooO'Or- 
O'O'O'O'O'OOOOOOOOOO'-'— ^r-r-»—*“f-(Mcgr>gcsJcvjCMCM 

oooooooooooooooooooooooooooooooo 



^ 00 < •- •• 

O II O M -) II O 

• • UJ t 

O “> O “> Q “> O ■ 

II O II O II O II c 

lO >0 I'- c 



in 



m 2 



• w 

r- r- c? in 

II II UJ II 

2C O' — 

+ 

o o o z: I 

o • o 3 m 

■o o <J m “j 

r- -o II II •> 

ns: z: s: 

0 in :3 0 3 ’- 

o “> i/> o m e: 



r- - o I 



195 
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ODIHENSIQN X( 72 ) , AL F A ( 72 ) , BETA ( 72 ) , G<\ MM A ( 72) , DELTA ( 72 ) ,UMB(72) , 02250 

1P(5,72),T(5,72) ,TAU(5,5) ,REC(5,5),R( 5,5) ,TX(5) ,TA( 10) ,CC (5) , 02260 
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RETURN 02530 

END 02540 



001 MENS ION X(72) , ALFA (72) , BETA! 72) , GAMMA (72) , DELTA (72) ,UMB( 72) , 02570 
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SUBROUTINE BVECTCR 02770 
OOIMENSION X(72) , ALFA (72), 6 ETA (72), GAMMA (72) , DELTA (72) ,UMB(72) , 02780 
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GO TO 1020 03040 

120 KI=22 03050 

NI=26 03060 

LI=21 03070 

1020 CONTINUE ,03080 
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1060 CONTINUE 0339C 

RETURN 03U00 

END 03410 
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GDI MENS ION X ( 72 ) , ALFA ( 72 ) , RET A ( 72 ) » GA MMA ( 72 ) , DEL T A ( 72 ) , UM B ( 72 ) , 0M94C 
1P(5,72),T(5,72) ,TAU(5,5) ,REC(5,5) ,R( 5,5) ,TX(5) ,TA(10),CC (5) , 04950 
2C(5,72),D(5,72) ,CK(5,72) ,SP(5,6) ,ST( 5,6) ,EXR(5,6), 04960 
3EA1 (72),UA(72 ) ,VA(72) ,WA(72) ,A(5C, 100) ,B(50) ,L{2C) ,X0{50 ) 04970 
OCOMMON X, ALFA, BETA, GAMMA, DELTA, P,T,TAU,REQ,R,ZZ,UMB, 04980 



ooooooooooooooooooooooo 
0'Oi-<Mr'OJtLnor-ooa'0*“<Mro^tD'Op>-ooo'0»- 
O'OOOOOOOOOO — '-r-r-'-*-'- — — f-tMCsi 

ooooooooooooooooooooooo 



OO'-OO — OOO-JOOOO. 



I O O O O 



>■ — r- + 

pn 

(<0 w r- W 

« d eg d 



in o 
o o 
eg eg 



II — II II II •— II II II II h- II II II II 
c\jcg^porr)i»i — inminin 

•— eg m > 1 - CM r 

— -'U.'-wwwu.wwww 

o m . o m 
eg ro • 



iJtiriwcgroj-in 



eg 



o m 
^ eg 
eg eg 



eg eg 



eg eg 



208 



A(6,5)=1.0 05220 

240 CONTINUE 05230 
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A(8,3) = A(9,3) 05550 
A(8,U) = A( 10, U) 05560 
A(8,5) = A( 1 1 ,5) 05570 
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An5,4J = USUM*(A( 10,4)-A(n,4) )/TU+A( 16,4) 05880 

A{15,5)=A( 16,5) 05890 
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An2,5) = A( 16,5) 06210 

06220 

CALL FIX (L,3) 06230 
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1 AC 26 , 3 ) 07180 

AC 25 , 4 )=USUM*(A( 20 , 4 )-A( 21 , 4 ) )/TU+A( 26 , 4 ) 07190 

A( 25 , 5 )=A( 26 , 5 ) 07200 



|0 



ooooooooooooooooooooooooooooooooooooo 

<\jfOj-vr»'Or~ooo'Oi“CMiviirunor-cocr'o<— (Njf<r»^i/>>or~oo 

oooooooooooooooooooooooooooooooo-ooooo 



cr 



o o o < < 

r<) M r<i I I 

3: r? X r- OJ 

»- t- 1 - « » 

► .. ► — ~ CC 00 

ro r- rr» o 00 •- t- 

•- « •- m w ^ 

^ ..•<t < — 

$ —o^»si<l— >— ro 

^o•— • •« 

— 2. IT) 

_J_I»JC'0'Z) D C\J 

r-t /1 </1 w 

3 : <t 

ino'iooo+ + -f 

X 13 

-J-J-JOOl- h- h- 

-s. •s. 

JJCOJ-CNJCsJ^ — — 

w W w <_ r- t\| rv) . 

» •- 

• o o 

^ ^ ^ ^ ^ IVJ (SJ C>J 

fOf^PO*— w W 

W<t < < 

-I -J -J -I -I I I I 

oo'ooo*— CM ro 

oocoocr O' O' 

siziacsis;'- ~ *- 

— 

n ooootO(/)io~ *- — 'j-vn 

I 33 XM<# * * - - 

u •.••■•►s; a; sioo 

3 3<MC\J 
3i :s33jcot/i (/)^w 

a w^www< X 3<< 

u 2: s: 5: s; 5; II ii ii ii ii 

►. 3 3 o 3 3 — — >- 

ooco«/}(/JV)<— <— cMCMro^w 

c 

- ^ -j ~J -J u a in ^ in Jt ^ 

L -J -J -J _J J CM CM CM CM CM CM CM 
-I <<<<<WWWWWWW 

i. UUUUU<<<<<<< 



► >o X • 

- w H- z; 

► _J -3 

- •. N- on 

► ^ • > 
i in CO ■»■ 
► w « X 
>3 0 1- 



CM J- CM 

W W w r- CM 



w w w < . < 

3 3 3 1 I 

O O O — CM 

O O O 00 00 



» 3 X * ' 

► • s; 

J < < 3 

> 3 3 <n 
- W w X 

•: s; s: II 
0 3 3 ^' 
n m m f- r 



CM ro m 






3 3 CM CM r 

< < w w ^ 

u u < < < 



3 •• 

H- ro 

CM 

» < 

o + 

► 3 



lO - 



3- m 



> w w w q: 



3*- 

•> SL ro jj- in >o 

< 3 CM CM C>i CM 

Z) I/) 

W 3 < < < < 

z: II II II II II 
3 — ^ — ^ — 

m r- CM ro J- m 

3 CM CM CM CM CM 
3 CM CM CM CM CM 

< W W w w W 

U < < < < < 



215 



oooooooooooooooooooooooooooooooooooooo 
O'— rjroijuior-coo'Oi— csjro-a-in'Of^oocror- csiroj-in'Ot'-oocT'Or- 
'00'0'0-o>0'0'0'0or^i^r^r-r^^-t^r~r-r^cocx3cocucoa5c3oooa3coo'cr'(7'cr'c7'a'c'a' 

r'-f-h-r-r'-r-r-r-r-r'-r-r^r-r-r-r-r'-r'-f'-h-t^r'-h-r'-h-r-r^h-i^r-f^r^r'-h-h-rw^.f*. 

OOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOO 



f 



ca o ki '- 

•^ -j 

< i<: (I •- 
-1 “3 II 
It in ^ 



•* in 

cu u ^ •> 

— o 

• 

< ^ II i- 
-J — II 
II H m Si 
z -J + -3 — m 



• • <t 

< II 
-3 “3 ~ 
II in 



< id 
•■3 “3 
II II 



ac. II II u I jc. II II u I ju II II u I 

_J -3 Z-3-I--3 O Z-3 + -3 — m Z-3-f O O 

< — in; 3 - “3 Oi— -3 I— f-J- “3 '-<MJ- 

Z 0 -I-+-I -00 **0 + + -l -00 ► 0 + + +00 •• -f + i-OC 
u o 3 <:i<:acoc 3 in'-*i<: 3 ^:*:'-'-incMi^ 3 (:a(:csj(Min m ic fo k 

o in“3“3“3intn + in"3“3~3inin + in“3T"3inin + < >“3"3int/ 

< II II II Z II II II Z II II II Z II II II II 

I— o<eouoo — o<o3ooo — o<o3uoo*-ic;<a3uc3C 

O O “3 ”3 -> O O < a -3 “3 -3 Q O < O -3 “3 “3 a Q < -3 -3 -3 "3 O C 



o 

o 

in 



o 



in 






in 



m o 

CM • 



cr — + 

m CM O' “3 -3 

ic: o ^ ^ 

o II (M II II 

«t no — 

-I 3 ^ “3 a < si -3 

m 

' $ 

a 



u 



o 



216 



ooooooooo 

0 ' 0 *-'(NirO; 3 ^tri*Oh- 

O'OOOOOOOO 

h-OOCOQOCDCOCOCDCD 

OOOOOOOOO 



ooooooooo 

ooo'o-— 

COOOCOCOCOCOOOCDCD 

ooooooooo 



oooooooo 

f^ooO'Or-rgroj' 

rgcjcMojoi 

cooooocooooooooo 

oooooooo 



ooooooooo 

in' 0 i^cD 0 'O<-rsiP 0 

(NjojCNjrsjcjm.'orofO 

cocococooococoooco 

ooooooooo 



•o tr» 



^ z z 



CC r- W ^ 



II < < < < < 

-J Z M II II II II V 
< ~ ~ ~ : 



II II II II II 



K in 



ro m • 

o - 

-« J- o 

CM m 



t m 



(T CM O r 

»- » in f 

W>0»- 

_J II «— o . 

+ Z II • -r II II • _j II II • -j II II • 

ZO~ZZO + ZZO + ZZOLU 
II II II — HO 

»-ir»in — »-inin— >h-inin-'minin— «z 



Jinf-o^^ — 0 "fn*— o 



c < < < < 



K h- 00 c 

Z CM CM C 

o — — ^ 

o < < « 



3 — -J O O 



U I U I — ^ r“ U I W I I U I U I a 

CMCMZ — roroZ»-^.:»Z*- — 

_| >0 M3 •• -J O >0 •• -J O >0 - h- O 

ro — Z Z— Z— ZZI-Q 

^ U.OO^U.OO^U.OO^U.00^0 IXJZ 

< MOO<Moa<Moo<MOo<o oc Lu 

o in o 



»- f- CM 



CM 



ro 



J- in 

O M> 



217 



ipwiiioa 



SUBROUTINE FIX(L,K) 08350 



oooooooooooooo 
'or-ooo' 0 »— cNi'viirinor^cocj' 

oooooooooooooooocoooooaocooo 

oooooooooooooo 



o oooooooooooooo 

•— csiro^in«or-a30'0'— CMfOjriD 
m in lA into inininm« 0 ' 0 > 0 ' 0 ' 0 <o 
00 oooooooooooooooooooooooooooo 
o oooooooooooooo 





► • — • o + 



z o r- : 

n ^ ^ 

— 2 : -I *-■ <— 

fcA eg II <-i 

Z I in~0'- 
LU ^ ro z m -J 

2 ; w. — ^ 



O r- Z 

>-« uj r- II — Lij 

— o — Z -J 3 

>j z in -J z z 
II I m II >-• ci: 
at o — H- o 
Z — Z Z I- Q 

— OU-O — OLUZ 
.J0»^0-J0QtUJ 



0000 



Z-Dr 
iO - 
O 3 »■ 



— 00 

► O' r- + 

o - -1 

I- ro o rj z 

— — r*- o r'- I 

I — . Z3 

-j y- Z3 r- a a. 

Z -J Z II I- II 

— — ro ro 
u. UL a. 3: c z 

tM K> o 



D 3 
z z 
— II 
u_ =? 



oino 00 omo 

fofo*- <M»n 'O'Or* 



o eg »© K ro in 
«o CO r«> r» h* N- 



218 



ooooooooooooooooooooooooooooo 
'Oh'Coo'Or-cMfOjtf^.inor-ooo'O — CMfna-unor-cDC'O'— CMfO^ 
•o -o -o -o r- r>- ^-^“^-^-^-oocococx^cocx^ooooc^coc^c^o^o'c^' 

OOOOCOOOOOCOOOCOOO CD 00 OOOOCOOOCOCOOOCOOOCOOOOOOOOOOOOOOOGO 

ooooooooOooooocooooooooooooooo 



o 

o 

ro 



00 — *o _J 

CO t-< 00 “> 



• — 

► CO -3 < h- < 

O CO • 3: CO 3: 

r- -j + •¥ 

— — -JSO — o zo 

•o ^rn_jiii 3 — o*-ro: 3 i- 

(VI (>-»— -3i-i(/)»-ro«— <— «/>*- 
I +13: I + 3t . 

_jcro-J=>tn II u-jo-j m a 

— II II — 3 w » 3 

u.^o-iu.ouoou.a-j(/)0 

-J O “3 I-I Q 26 o I-I o “3 3: o 



ro -3 
O' — 

•“ <; 

t\j 5: 

o + 

• z 

CM 3 
•— O' to , 

+ -' 3: 

_l 3 II o 3 CO I 

z z z M h- z O' : 

II — 3 ro II : 
-J U. (O 2 O 3 a t 

3 •-« 3: h» t5 3 o ; 






► _J 

O (O 3 
I- O CM 

CM ro »— < 

•— ► * 3C 

► O (M 

O CM CM + 

CM *— r- 

• z 



• -J CM 3 

O 3 CM to 

CM I- r- 3. 




-J -J <-« _l z 

w (I ^ 3 

U. _J U. LL CO 

3 3S 



o cMf'-'Oco 

N- 00 op r— 00 CO 00 



O <M fO lA 

O' CJ> O' O' 



O »- O CM 

r- CM O CM 



219 



TWl = 1 08950 
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SUBROUTINE LINEAR ( A , B , MRCW , NCCL , XO , YC ) 09290 
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DO 1340 1= 1 , 8 OOnO 

TERR(I) = 0.0 C0120 

IOF[X( 1+8) = 0 C0130 

1340 lOFIXU) = INFIX(I) 00140 
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1043 IF (lA) 1050, 1048, 1C5C 01600 

1048 IF ( Y(I) ) 1050, 105C, 1045 01610 

1044 IF (lA) 1045, 1046, 1045 01620 

1045 IF ( YI - AA ) 1C50, 1050, 1047 01630 

1046 lA = 1 01640 
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Y(IRJ = -1. 01960 

LL = 0 01970 

903 00 904 L = IR, M2, M 01980 

IF ( E(L) ) 905, 914, 905 01990 

914 LL = LL + M 02000 
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LL = (JM - n * ME 02670 

301 DO 303 MM = 1, M 02680 

LL = LL + 1 02690 

IF ( A( LL ))304, 303, 3C4 02700 

304 DT = DT + P( MM J * A ( LL ) 02710 
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